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Let  7Yinf  be  the  family  of  functions  each  of  which  is  an  infimum  of  arbitrary 
collection  of  positive  harmonic  functions  on  a  bounded  domain  D  in  Rn.  Hild  is  a 
convex  cone  of  positive  functions  closed  in  the  topology  of  pointwise  convergence  and 
closed  under  arbitrary  infima.  Every  such  function  is  excessive  for  Brownian  motion 
in  D  and  satisfies  Harnack  inequality. 

The  most  important  functions  in  the  class  are  those  defined  by 

sa  =  mf{h  :   h  >  0  harmonic,  h(a)  =  1},    a  €  D, 

which  give  the  best  estimate  for  the  constant  in  Harnack  inequality.  Those  functions 
are  explicitly  computed  for  the  unit  ball  and  the  halfspace.  It  is  shown  for  Lipschitz 
domains,  that  the  knowledge  of  functions  sa  is  enough  to  recover  all  positive  harmonic 

functions.  More  precisely, 

r        sa(x) 

lim  — 7— ^-  =  A  (x,  z) 

«-»*  sa(x0) 

for  every  z  £  3D,  where  K  is  the  kernel  function  for  D  at  x0. 

Several  examples  of  functions  in  H™1  are  given;  for  instance,  the  expected  exit 
time  from  D  for  the  Brownian  motion  is  a  member  of  the  family.  A  necessary  condi- 
tion satisfied  by  the  Riesz  measure  of  any  function  in  Hini  is  obtained. 


IV 


Results  are  proved  in  a  more  general  setting  which  makes  it  possible  to  simulta- 
neously treat  both  harmonic  functions  and  positive  solutions  of  the  equation  Lu  =  0, 
where  L  includes  uniformly  elliptic  operators  in  divergence  form  or  Schrodinger  op- 
erators. Positive  harmonic  functions  for  a  certain  class  of  diffusions  are  also  included 
(with  more  restrictive  results). 

It  is  shown  that  positive  harmonic  functions  determine  Brownian  motion  up  to 
time  change.  Moreover,  any  diffusion  in  a  Lipschitz  domain  D  whose  killing  distri- 
butions are  equal  to  the  exit  distributions  of  Brownian  motion  from  D  is  necessarily 
time  changed  Brownian  motion. 


CHAPTER  1 
INTRODUCTION 


Positive  harmonic  functions  on  a  bounded  domain  D  in  Rn  enjoy  two  impor- 
tant and  well-known  properties:  They  satisfy  Harnack  inequality  and  allow  integral 
representation  via  a  kernel  function. 

Harnack  inequality  is  usually  expressed  in  two  forms.  The  first  one  is  more  local 
in  nature  and  states  that  if  a  ball  B(xo,r)  is  contained  in  D,  then  for  any  positive 
harmonic  function  h  on  D,  and  any  x  G  B(x0,r), 

(r  +  |x  —  a;0|jn  (r  —  |x  —  x0|Jn 

The  second  form  follows  from  the  inequality  above:  For  any  compact  subset  K 
of  D  there  is  a  constant  c  =  c^  such  that  for  any  positive  harmonic  function  h  on  D 

Tr\<c    for  all   x,y  €  K  (1.2) 

%) 

(see,  e.g.,  [14,  p.  31]). 

Harnack  inequality  has  numerous  applications.  We  concentrate  on  the  following 
two:  (i)  the  infimum  of  any  family  of  positive  harmonic  functions  satisfies  both  (1.1) 
and  (1.2),  and  (ii)  any  such  infimum  is  continuous.  The  second  property  follows  from 
the  fact  that  every  function  satisfying  (1.1)  is  automatically  continuous  at  x. 

Hence,  one  can  form  the  family 

Ji      =  {u  :  u  =  mfha,  ha    positive  and  harmonic  in  D} 

consisting  of  continuous  functions  satisfying  Harnack  inequality.    Every  function  in 

Hivi  is  trivially  superaveraging,  so  H"**  consists  of  continuous  positive  superharmonic 
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functions  on  D  which  satisfy  Harnack  inequality.  Another  way  of  expressing  this  is 
that  functions  in  7iild  are  excessive  for  the  Brownian  motion  in  D  and  satisfy  Harnack 
inequality. 

In  Chapter  2  we  prove  the  following  result  about  HinI: 

THEOREM  2.1  7ilnt  is  a  convex  cone  stable  for  arbitrary  infima  and  closed  in  the 
topology  of  pointwise  convergence. 

Now  we  recall  the  second  essential  property  of  positive  harmonic  functions, 
namely  the  integral  representation. 

Following  [11]  let  D  be  a  domain  in  Rn  with  the  Green  function  Go  and  v  a 
measure  with  compact  support  Suppi/  C  D.  The  function 

GDv{y) 

is  called  the  Martin  function  (or  the  Martin  kernel  )  based  on  v.  There  is  a  unique 
metrizable  compactification  Dm  of  D  such  that  each  Martin  function  Ku  has  a  contin- 
uous extension  (denoted  also  by  K„)  to  D  x  (Dm  \  Supp^),  and  K„(-,y\)  —  Ku(-,y2) 
if  and  only  if  j/i  =  y2.  The  boundary  &mD  =  Dm  \  D  is  called  the  Martin  boundary. 
Let  dMD  denote  the  set  of  minimal  points  z  G  OmD.  Then  for  each  positive  harmonic 
function  h  on  D  there  is  a  unique  measure  \i  on  dMD  such  that 

h=  J        K„(;z)n(dz).  (1.3) 

Jd"MD 

If  the  measure  v  is  the  point  mass  at  x0  €  D,  then  we  say  that  the  Martin 
kernel  is  based  at  x0.  In  this  case  the  continuous  extension  of  the  Martin  kernel  to 
8mD  means  that  as  y  — »  z  in  the  Martin  topology,  where  y  €  D,  z  G  6mD,  we  have 


For  the  unit  ball  D  =  5(0,1)  the  Martin  boundary  is  simply  the  Euclidean 
boundary  of  D,  and  the  Martin  kernel  is  the  Poisson  kernel 

1-lzl2 


P(x,z)  = 


x 


Inequality  (1.1)  for  the  unit  ball 


MO)  <  h(x)  <  /      |g|.2    h(0) 


(1  +  |*|)»  w"  v  '-  (l-\x\)" 
reveals  the  connection  with  the  Poisson  kernel.  In  fact,  it  is  derived  from  the  repre- 
sentation formula.  Moreover,  this  inequality  is  sharp  in  the  sense  that  the  bounds  are 
attained  for  h  =  P(-,z)  where  z  =  x/\x\  or  — ar/|x|.  In  a  sense,  Harnack  inequality  is 
tailor-made  for  balls.  It  is  not  optimal  for  other  domains. 

To  get  the  optimal  lower  bounds  for  an  arbitrary  domain  D,  it  is  natural  to 
introduce  the  function 

k(x,y)  =  mi{h(x)  :  h  positive  harmonic  in    D,  h(y)  =  1},     y  E  D.  (1.5) 

From  Theorem  1  this  function  is  continuous  in  x  (in  fact,  it  is  jointly  continuous;  see 
Proposition  2.1)  and  gives  the  greatest  lower  bound  for  a  positive  harmonic  function 
which  is  1  at  y.  As  in  the  case  of  the  unit  ball,  this  lower  bound  is  attained  by  a 
Martin  kernel:  There  exists  z  —  z(x,y)  in  the  minimal  Martin  boundary  d°MD  such 

that 

u       ^       K(x'z) 
K{y,z) 

where  K  is  based  at  some  point  x0  in  D  (see  Proposition  2.4).  Therefore,  the  Martin 
kernel  K(x,  z)  completely  determines  function  k. 

For  Lipschitz  domains  the  converse  is  also  true:  Function  k  is  sufficient  to 
recover  the  Martin  kernel.  More  precisely,  let  D  be  a  bounded  Lipschitz  domain. 
Then  the  Martin  boundary  dMD  of  D  is  exactly  the  Euclidean  boundary  (see  [15]), 


and  all  boundary  points  are  minimal.  Let  K  be  the  Martin  kernel  based  at  x0  G  D. 
In  Theorem  2.2  we  prove  that  if  x  G  D,  and  z  €  3D,  then 

lim^Mr  =  #(*,*).  (1.6) 

This  result  can  be  regarded  as  an  analogue  of  (1.3).  The  reason  we  restrict  ourselves 
to  Lipschitz  domains  is  that  our  proof  relies  on  the  fact  that  if  Zi,Z2  £  dD,  z\  ^  z25 
x  €  D,  z  €  dD,  then 

lim       K(x,z)  =  0.  (1.7) 

This  is  not  true  in  general:  There  are  examples  of  domains  such  that  zi,  z2  are 
distinct  minimal  boundary  points  on  the  Martin  boundary,  but  limA'(x,z)  =  +oo 
as  (x,z)  — >  (zi,Z2)  (see  [19]).  We  note  that  (1.7)  is  also  true  for  so-called  "NTA" 
domains  introduced  in  [16].  That  the  limit  in  (1.7)  is  zero  follows  from  Lemma  5.2 
in  [16,  p.103]. 

In  Sections  4.1  and  4.2  we  obtain  explicit  formulae  for  k  in  case  D  is  the  unit 
ball  or  the  upper  halfspace.  We  note  that  the  similar  formulae  were  also  obtained  in 
[17].  Computations  rely  on  the  Mobius  and  Kelvin  transform.  Those  formulae  show 
that  k  is  not  symmetric  if  dimension  n  >  3.  For  n  =  2,  k  is  symmetric,  which  is  also 
true  for  some  other  domains  in  R2. 

As  we  have  mentioned  above,  7iini  is  a  convex  cone  of  positive  functions.  It 
is  natural  to  ask  whether  7{inf  is  a  potential  cone.  Although  general  wisdom  tells 
us  the  answer  is  no,  it  does  not  seem  that  potential-theoretic  arguments  can  justify 
that,  simply  because  in  case  n  =  1  the  answer  is  yes.  In  Section  4.4  we  give  a  simple 
example  for  n  =  2  that  confirms  that  no  potential  theory  is  incorporated  in  Hird. 

An  intriguing  question  to  us  is,  given  a  continuous  excessive  function  on  D, 
how  can  we  recognize  it  as  an  infimum  of  positive  harmonic  functions?  In  Section 
3.3  we  give  a  partial  answer:    If  u  =  Gd(*  is  a  potential  in  ?iinf,  then  //  cannot 


have  a  compact  support  (provided  that  D  "has  no  holes").  This  is,  of  course,  very 
far  from  being  sufficient.  We  do  not  know  if  any  condition  on  a  measure  fi  can  be 
given  to  guarantee  that  £?£>//  is  in  Hlni.  The  next  example  should  provide  some 
insight  into  the  difficulties.  Let  D  be  the  unit  disc  in  R2,  z  =  (1,0)  €  dD,  and 
u  =  P(-,z)  A  1.  We  show  in  Section  4.5  that  u  is  extremal  in  H1"1.  It  is  easy  to  see 
that  the  Riesz  measure  of  u  is  concentrated  on  the  boundary  S  of  the  disc  B(xQ,  1/2) 
where  xq  —  (1/2,0).  Let  A,B  form  any  partition  of  S  such  that  neither  A  nor  B  is 
compact.  Then  u  =  G&\i\a  +  GdH\b-,  and,  since  u  is  extremal  in  7iini,  at  least  one  of 
the  potentials  GdH\a,  Gd\i\b  is  not  in  7iinf. 

In  view  of  the  above,  it  is  not  surprising  that  not  many  instances  of  functions 
in  7iinf  are  available.  We  give  several  examples  in  Section  4.3,  the  most  interesting 
being  the  expected  exit  time  of  the  Brownian  motion  from  D. 

The  organization  of  the  thesis  is  directed  by  an  attempt  to  include  more  general 
situations.  Topological  properties  of  the  cone  TV**  depend  exclusively  on  the  fact  that 
positive  harmonic  functions  on  a  domain  form  a  closed  convex  cone  with  compact 
basis.  Therefore,  in  Chapter  2,  we  study  closed  convex  cone  H  with  compact  basis  of 
positive  continuous  functions  on  a  locally  compact  space.  We  introduce  the  function 
k  as  in  (1.5)  and  the  cone  H,nI  of  all  infima  of  function  from  H.  We  note  that  a 
similar  approach  can  be  found  in  [7]. 

In  Section  2.2  we  assume  the  existence  of  a  kernel  function  K.  We  require  K 
to  have  properties  of  the  Martin  kernel  on  a  Lipschitz  domain  and  postulate  them  in 
hypotheses  (Hi)-(HZ).  In  this  context  we  prove  that  (1.6)  holds. 

In  Chapter  3  we  give  examples  of  cones  with  compact  basis.  Instead  of  positive 
harmonic  functions  for  the  Laplacian,  one  can  consider  positive  solutions  to  Au  =  0, 
where  A  is  a  uniformly  elliptic  operator  in  divergence  form.  By  using  results  in  [4] 
which  parallel  those  obtained  in  [15],  we  prove  that  (1.6)  is  valid,  where  now  K  is 


the  kernel  function  for  (A,D),  D  Lipschitz,  and  k  denned  by  (1.5)  for  the  cone  of 
positive  harmonic  functions  for  A. 

In  Section  3.2  the  same  result  is  obtained  for  positive  solutions  of  Schrodinger 
equation  —  Au  +  qu  =  0,  for  q  in  Kato  class  with  finite  gauge  and  D  Lipschitz.  Here 
we  use  results  from  [5]  and  [9].  We  note  that  equivalent  statements  can  be  obtained 
for  positive  solutions  of  the  equation  ^Au  +  b  •  Vu  =  0  by  using  results  from  [10]. 
The  drift  b  must  satisfy  conditions  given  in  that  paper. 

In  Section  3.3  we  use  probabilistic  methods  and  show  that  positive  harmonic 
functions  for  certain  diffusions  also  form  a  convex  cone  with  compact  basis.  We  rely 
on  results  from  [21]  and  [22].  Of  course,  in  this  setting  the  kernel  function  satisfying 
our  hypotheses  (Hi)-(H3)  is  not  available. 

In  Chapter  4  we  return  to  classical  harmonic  functions  and  obtain  some  formulae 
and  examples  described  above.  In  addition,  we  obtain  some  (but  not  all)  extremal 
elements  in  Ti*"*  on  the  unit  ball. 

In  Chapter  5  we  investigate  a  diffusion  on  a  bounded  domain  D  which  has  the 
same  class  of  harmonic  functions  as  the  Brownian  motion  in  D.  We  show  that  such 
a  diffusion  is  necessarily  a  time  change  of  the  Brownian  motion.  Moreover,  if  the 
diffusion  dies  as  it  approaches  the  boundary  3D  with  the  same  distribution  as  the 
Brownian  motion,  it  is  again  a  time  change. 


CHAPTER  2 
CONVEX  CONES  OF  POSITIVE  CONTINUOUS  FUNCTION 


2.1     Topological  Properties 

Let  D  be  a  locally  compact  topological  space  with  countable  basis.  Then  D 
is  metrizable  and  let  d  denote  a  metric  compatible  with  the  topology  on  D.  Let 
{Kn}^=1  be  an  increasing  sequence  of  compact  sets  covering  D  such  that  Kn  lies  in 
the  interior  of  Kn+i  for  each  n  G  N. 

Let  C(D)  denote  the  vector  space  of  all  real- valued  continuous  functions  on  D, 
topologized  by  the  separating  family  of  seminorms 

pn(f)  =  aup{|/(s)| :  x  6  Kn). 

The  topology  on  C(D)'\s  metrizable,  and  C(D)  thus  becomes  a  Frechet  space. 
Convergence  in  this  topology  is  the  uniform  convergence  on  compact  subsets  of  D. 
When  we  refer  to  any  topological  property  of  C(D) ,  we  always  refer  to  this  topology. 

A  subset  B  of  C(D)  is  bounded  if  and  only  if  for  each  compact  subset  K  of  D, 
there  exists  a  constant  CK  such  that  |/(a:)|  <  CK  for  each  f  £  B  and  each  x  G  K. 

Let  A  C  D  and  U  C  C(D) .  U  is  said  to  be  uniformly  equicontinuous  on  A  if 
for  every  e  >  0  there  exists  S  >  0  such  that  \u(x)  —  u(y)\  <  e  whenever  d(x,y)  <  S, 
x,y  G  A,  u  G  U.  A  family  U  C  C(D)  is  locally  uniformly  equicontinuous  if  it  is 
uniformly  equicontinuous  on  each  compact  subset  K  of  D. 

We  shall  freely  use  Arzela-Ascoli  theorem:  A  family  U  C  C{D)  is  compact  if 
and  only  if  it  is  closed,  bounded,  and  locally  uniformly  equicontinuous. 
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By  C+(D)  we  denote  the  family  of  strictly  positive  functions  in  D.  Let 
H'  C  C+(D)   and  let  W  =  H'U  {0},  where  0  denotes  the  function  which  is  identically 
zero. 

We  assume  that  Ti  is  a  convex  cone  of  functions,  closed  in  C{D) ,  having  a 
compact  basis.  The  last  condition  means  that  there  exists  a  hyperplane  C  in  C(D)  such 
that  the  family  Q  =  C  C\  Ti'  is  compact  in  C(D)  and  generates  H: 

H-  {Xu:  A>0,ue£}. 

The  hyperplane  C  =  {/  €  C(D)    :  <  <f>,f  >=  a},  where  <f>  is  a  continuous  linear 
functional  and  a  €  R. 

Having  established  the  setting,  we  proceed  by  defining  the  family  of  all  possible 
infima  of  functions  in  H: 

HinS  =  {u  :u  =  inf  ua,ua  <E  H}.  (2.1) 

The  main  purpose  of  this  section  is  to  show  that  the  family  7iinfconsists  of 
strictly  positive  continuous  functions  (including  zero)  and  that  it  is  closed  in  the 
topology  of  pointwise  convergence.  The  basic  tool  toward  this  goal  is  the  function 
k-.DxD^R  defined  by 

*(*.*)  =  $,g.  (2.2) 

This  function  gives  the  best  estimate  for  u(x)  when  u(y)  has  some  prescribed 
value.  This  will  be  made  more  precise  later.  We  emphasize  that,  in  general,  k  is  not 
symmetric  in  x  and  y. 

Note  that 

•   r   u(x)       .   ,  u(x) 

3e*)m2b&  (23) 

hence  k  could  be  computed  using  only  functions  in  the  basis  Q. 
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Further,  for  x,  y,  z  in  Z), 

*(„,)*<„„)  -  *r,g  -%;M  <  ^M  .  k(x,y).     (,4) 

By  taking  x  =  y  and  using  the  trivial  fact  that  k(x,x)  =  1,  we  get 

k(x,z)k(z,x)  <  1,     x,z(zD.  (2-5) 

LEMMA  2.1  For  all  x,y  in  D,  k(x,y)  >  0. 

u{  x ) 
PROOF:    By  (2.3),  k(x,y)  =  inf  Hence,  there  is  a  sequence  {un}  C  £/  such 

uzQ  u{y) 

that  un(x)/un(y)  — ►  k(x,y).  By  compactness  of  Q,  there  exist  a  subsequence  {un)} 
and  u  E  Q  such  that  un   — »  u  in  C(Z)) .  Therefore, 

k(x,y)  =  hm  — ^  =  -^  >  0, 

since  u  >  O.d 

Since  all  functions  in  W  are  strictly  positive,  one  can  consider  the  family  Q~x  = 

{1/u  :  u  €  Q}  which  is  also  contained  in  C+(D)  .  This  family  is  compact  in  C(D) . 

Indeed,  given  a  sequence  {l/un}  in  G~x,  compactness  of  Q  yields  a  subsequence  {un  } 

and  u  €  Q  such  that  unj  — ►  u  in  C(D) .  Hence,  l/unj  -»  l/u  pointwise,  and  uniform 

convergence  follows  from 

|   1         1,  _  W-unj\ 

Un,  U  \unj\\u\ 

and  boundedness  of  the  subsequence  on  compacts. 

Let  K  be  a  compact  subset  of  D.  Since  Q  is  bounded,  there  is  a  constant 
C  =  Ck  such  that  u{x)  <  C  for  each  u  G  Q  and  each  x  6  K.  Similarly,  there  is 
c  =  cK  such  that  l/u(x)  <  1/c  for  each  l/u  e  £_1  and  for  each  x  e  K.  Hence, 

c  <  u(z)  <  C,       u€G,  x€K.  (2.6) 
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LEMMA  2.2   For  each  y  6  A",  the  function  x  \~*  k(x,y)  is  uniformly  continuous  on 
K. 

PROOF:  Let  e  >  0.  By  the  local  uniform  equicontinuity  of  Q,  there  exists  8  >  0  such 
that 

\u(x)  —  u(x')\  <  ce,    x,x'  £  K,  d(x,x')  <  8,  u  £  Q, 

where  c  is  the  constant  from  (2.6).    This  can  be  written  as  u(x)  <  u(x')  +  ce  and 
u(x')  <  u(x)  +  ce.  After  dividing  the  first  inequality  by  u(y)  and  using  c/u(y)  <  1, 

we  get 

u(x)       u(x') 
—^-L  <  — — -  +  e. 
u(y)        u(y) 

By  taking  the  infimum  over  all  u  in  Q,  it  follows  k(x,y)  <  k(x',y)  +  e.    Similarly, 

k(x',y)  <  k(x,y)  +  e,  i.e., 

|^(x,  y)  -  k(x',  y)\<e,    x,  x  £  K,  d(x, x')  <  8.  (2.7) 


LEMMA  2.3  For  each  x  £  D,  the  function  y  i— >  k(x,y)  is  continuous  at  x. 

PROOF:    Without  loss  of  generality,  let  us  assume  that  x  £  IntA',  where  K  is  the 
compact  set  as  above.  Let  e  >  0  and  choose  8  >  0  such  that    \u(x)  —  u(y)\  <  ce, 
y  £  K,  d(x,y)  <  8,  u  £  Q.  After  dividing  by  u(y),  it  follows  that  1  -e  <  u(x)/u(y)  < 
1  +  e,  y  £  K,  d(x,y)  <  8,   u  £  Q.    By  taking  the  infimum  over  all  u  £  Q,  we  get 
1  —  e  <  k(x,y)  <  1  +  e.  Since  k(x,x)  =  1,  the  above  reads 

\k(x,y)  -  k(x,x)\  <  e,    y  £  K,  d(x,y)  <  8.  (2.8) 


□ 
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Now  we  are  ready  to  establish  the  joint  continuity  of  k. 

Let  (x,y)  €  D  x  D  and  let  {(xn,yn)}^=1  be  a  sequence  in  D  x  D  converging  to 
(x,y).  Then, 


From  (2.4), 


Therefore, 


Hxn,yn)  __  Hxn,yn)k(xn,y) 
k(x,y)  k(xn,y)    k(x,y)  ' 


it         \  ^  k{xn, yn)  1 

k{y,yn)  <  77 r  < 


Kxn,y)  "  Hyn,y)' 


Hv,v.)kp^  <  !p*£  <    i    MM  (2.9) 

fc(x,y)  «(ar,y)         «(y»,  t/)  *(*>y) 

As  n  — >  oo,  k(y,yn)  — >  k(y,y)  by  Lemma  2.3,  A;(xn,j/)  — >  k(x,y)  by  Lemma  2.2.  Let 
n  — >  oo  in  (2.9);  it  follows 

1  <  hminf  — p; r1-  <  hmsup  -7- -^  <  1, 

«(x,y)  n         «(x,y) 

which  shows  that  lim  k(xn,  yn)  =  k(x,y).  Thus,  we  have  proved  the  following. 

PROPOSITION  2.1    The  function  k  :  D  x  D  ->  R  de/med  6y  (2.2)  is  s*n'c%  positive 
and  (jointly)  continuous. 

A  family  of  nonnegative  functions  on  D  is  usually  said  to  satisfy  Harnack  prin- 
ciple if,  for  every  compact  subset  K  of  D,  there  is  a  constant  C  =  CK  such  that 

u(x)  <  Cu(y)  (2.10) 

for  all  x  and  y  in  K,  and  all  u  in  the  family.  The  proposition  above  trivially  implies 
the  following 

COROLLARY  2.1  Let  H  be  a  closed  convex  cone  as  above.   Then  H  satisfies  Harnack 
principle. 
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PROOF:  If  K  is  compact,  then  the  function  k  attains  its  strict  positive  minimum, 
say  m,  on  K  x  K.  Then  u(x)/u(y)  >  k(x,y)  >  m  for  all  x,y  €  K  and  any  u  £  H.n 
Now  we  slightly  change  our  point  of  view.  Instead  of  getting  the  function  k 
from  a  given  cone  7i,  let  us  assume  that  we  are  given  a  (jointly)  continuous  function 
k  :  D  x  D  — »  (0,oo)  which  is  identically  1  on  the  diagonal,  i.e.,  k(x,x)  =  1  for  each 
x  €  D.  Using  k,  we  define  the  family  S  of  all  nonnegative  functions  u  :  D  — ►  R 
satisfying 

u(x)>  k(x,y)u(y)  (2.11) 

for  all  x,y  G  D. 

The  following  lemma  contains  some  properties  of  the  family  S. 

LEMMA  2.4  (i)  S  is  a  convex  cone  of  continuous,  strictly  positive  functions  (unless 
identically  zero). 

(ii)  S  is  closed  in  the  topology  of  pointwise  convergence  (and,  hence,  in  C(D) ). 

(hi)  S  is  stable  under  arbitrary  infima  and  suprema. 

(iv)  IfTcS  is  bounded  at  some  point  x  £  D,  then  T  is  bounded  in  C(D) . 

(v)  IfTdS  is  bounded  in  C(D) ,  then  it  is  locally  uniformly  equicontinuous. 

(vi)  Let  x  G  D  and  Sx  =  {u  e  S  :  u(x)  =  1}.   Then  Sx  is  a  compact  basis  for  S. 

PROOF:  (i)  If  xn  -►  x,  then  liininf  u(xn)  >  liminf  k(xn,x)u(x)  =  u(x) 

>  limsupnfc(x,a;n)u(xn)  =  limsupnt<(xn),  which  gives  continuity.  The  rest  of  (i),  as 

well  as  (ii)  and  (iii)  is  trivial. 
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(iv)  Let  C  be  a  constant  such  that  u{x)  <  C  for  each  u  €  T.  Let  K  be  a  compact 
set  in  D  and  a  =  inf  k(x,  y)  >  0.  Then,  for  every  y  G  K, 

C  >  u(x)  >  k(x,y)u(y)  >  au(y), 

and  so  u(y)  <  C/a,  which  proves  that  T  is  bounded  in  C(D) . 

(v)  Let  C  be  a  constant  such  that  u(x)  <  C  for  each  u  €  T  and  each  x  G  K.  For 

e  >  0,  there  is  6  >  0  such  that,  for  all  (xi,yi)  and  (x2,t/2)  in  K  x  /i'  with 

di{(xi,yi),(x2,yi))  =  d(x1,x2)  +  d(yuy2)  <  S,  we  have  |fc(x1;j/i)  -  k(x2,y2)\  <  e/C. 

Hence,  for  x,y  £  K  with  d(x,y)  <  8,     \k(x,y)  —  1|  =  |A;(x,t/)  —  A:(a:,x)|  <  e/C,  and 

|Ar(j/,x)-l|<e/C. 

For  mGT, 

u(i/)  -  «(x)  <  u(y)  -  k(x,y)u(y)  =  (1  -  k(x,y))u(y), 

and 

u(x)  -  u(j/)  <  u(x)  -  fc(y,  x)u(y)  =  (1  -  fc(y,  x))«(x). 

If  ti(y)  >  u(x),  then  |u(y)  —  u(x)|  <  |1  —  k(x,y)\  u(y)  <  C\l  —  k(x,y)\.  Similarly  if 

u(x)  >  u(y),  then  \u(y)  —  u(x)\  <  C\\  —  k(y, x)\. 

Hence,  \u(y)  —  u(x)\  <  Cmax{|l  —  k(x,y)\,  |1  —  k(y,x)\}  <  e,  whenever  x,y  G  K, 

d(x,  y)  <  8,  and  mGT. 

(vi)  By  definition,  Sx  is  bounded  at  x.  Hence  (iv)  implies  boundedness  in  C(D) , 

which,  by  (v),  yields  local  uniform  equicontinuity.  Trivially,  Sx  is  closed,  so 

Arzela-Ascoli  theorem  gives  compactness.  That  Sx  is  a  basis  for  S  is  obvious.  □ 

To  recapitulate:  Starting  from  a  strictly  positive,  continuous  function  k  satisfy- 
ing k(x,x)  =  1,  we  have  obtained  a  closed  convex  cone  with  compact  basis  of  strictly 
positive  (unless  zero)  continuous  functions.  Some  care  must  be  exercised,  because  for 
some  functions  k  we  get  trivial  cones.  For  example,  if  k  is  symmetric  and  k(x,  y)  >  1 
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for  some  x,y  6  D,  the  only  function  which  satisfies  (2.11)  is  zero.  If  k  is  identically  1, 
we  obtain  only  nonnegative  constants.  On  the  other  hand,  if  k  also  satisfies  condition 
(2.4),  then  for  every  y  6  D,  the  function  x  i-»  k(x,y)  is  in  S.  A  non-trivial  example 
of  such  a  function  is  k(x,y)  =  exp(—d(x,y)). 

Let  T  be  a  closed  convex  cone  contained  in  S.  Then  Tx  =  {u  €  T  :  u(x)  =  1}  is 
a  compact  basis  for  T.  Besides  being  closed  in  the  topology  of  uniform  convergence 
on  compacts,  T  is  also  closed  in  the  pointwise  topology.  Indeed,  if  {un}^Li  is  a 
sequence  in  T  such  that  un(x)  — >  u(x)  for  every  i^D,  then  {un(x0)}^=l  is  a  bounded 
sequence,  x0  being  any  point  in  D.  By  (iv)  and  (v)  of  Lemma  2.4,  {un}^=1  is  bounded 
in  C(D)  and  locally  uniformly  equicontinuous,  hence  relatively  compact.  Therefore, 
there  exist  a  subsequence  {uUj}  and  a  function  v  €  C(D)  such  that  un  — >  v  in  C(D) . 
Since  T  is  closed,  v  €  T.  But  unj  — >  u  pointwise,  so  u  =  v  £  T. 

Let  us  now  go  back  to  the  original  setting:  7i  is  a  closed  convex  cone  of  strictly 
positive  (unless  zero)  continuous  functions,  having  a  compact  basis  Q.  The  function 
k  :  D  x  D  — >  (0,oo)  is  defined  as  in  (2.2).  Using  this  k,  we  form  the  cone  S  of 
functions  satisfying  (2.11).  Trivally,  u(x)  >  k(x,y)u(y)  for  each  u  G  H,  and  so 
HCS. 

Let  Ji  =  {uiAu2A.../\un:uj  eH,  j  =  l,2,...,n,  n  e  N}  and  J  =C1  Jx  be 
the  closure  of  2^  in  C(£>) .  Both  I1  and  J  are  contained  in  S.  J  is  a  closed  convex 
cone,  stable  under  finite  minima.  Convexity  follows  from  the  equality 

(tii  A  ...  Aun)  +  (vx  A  ...  Avm)  =  min{(w,-  +  Vj)  :  i  =  l,...,n,  j  —  l,...,m}. 

Obviously,  J  is  the  smallest  closed  convex  cone  stable  under  finite  minima 
containing  H.  The  family  lXo  =  {u  G  J  :  u(ar0)  =  1},  :r0  G  £>,  is  compact  and  J 
is  closed  for  pointwise  limits.  Furthermore,  I  is  stable  for  countable  infima.  Indeed, 
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if  {«n}£Li  is  a  sequence  in  J  and  u  =  inf  un,  let  vn  =  Ui  A  . . .  A  un  G  X.    Since 
u  =  lim„  vn,  u  is  in  X. 

Now  we  recall  that  the  family  of  all  possible  infima  of  functions  in  H  was  defined 
by  (2.1)  and  was  denoted  by  Winf.  Since  each  function  in  7i  satisfies  (2.11),  the  same 
is  true  for  functions  in  H,nI .  Hence  7i,nf  C  S  and,  in  particular,  each  function  in  Hini 
is  continuous. 

LEMMA  2.5  Each  function  in  HinS  is  an  infimum  of  a  countable  family  of  functions 
in  7i. 

PROOF:    Let  u  =    inf  ua,   ua  G  "H.    By  Choquet's  lemma  (e.g., [14,  p. 34]),  there 

exists  a  countable  set  Ao  C  A  with  the  property  that  if  g  is  a  lower  semicontinuous 

function  on  D  such  that  g  <    inf  ua,  then  g  <  u.  Take  g  =    inf   ua  and  note  that 

azAs  <*e«4o 

g  is  continuous.  Then    inf   ua  <  u,  and  since  the  reverse  inequality  trivially  holds, 

aeAo 

u  =    inf   ua.u 

aeAo 

Let  us  temporarily  introduce  the  following  notation: 

J  =  {u  :  D  ^  R.u  =lim«n,  m  <  u2  <  . . . ,  un  G  Hiai}. 

n  J 

Each  function  in  J  satisfies  (2.11),  so  J  C  S. 

In  the  following  three  lemmas,  we  show  that  all  three  families,  Hini,  J,  and  I 
are  equal. 

Lemma  2.6  J  =  1 

PROOF:  By  Lemma  2.5,  each  function  in  Hitd  is  a  countable  infimum  of  functions 
in  H,  hence  7illd  C  X.  Since  J  is  closed  for  pointwise  convergence,  it  follows  that 
J  CX. 
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Conversely,  let  u  (EX.  Then  u  =  limun,  where  un  6  J:  C  Hini.  Let  Vk  =  inf  un, 

n  n>k 

for  igN.    Then  Vk  €  X,  f *  £  "Winf,  and  {u,t}  is  an  increasing  sequence.    Hence, 


i>  =  \\mvk  €  »7.  But 
k 


so  u  €  J.  D 


u  =  sup  Vk  =  lim  inf  un  =  u, 
k  n 


LEMMA  2.7  Let  u  e  Hild  and  x  £  D.  Then  there  exists  v  e  H  such  that  v  >  u  in  D 
and  v(x)  =  u(x). 

PROOF:  Let  u  =  inf  uQ.  There  is  a  sequence  {un}  C  {ua}  such  that  un{x)  [  u(x). 
Thus  {un}  is  bounded  at  x,  and,  by  (iv)  and  (v)  of  Lemma  2.4,  it  is  relatively  compact. 
Therefore,  there  is  a  subsequence  {uU}}  and  v  £  H  such  that  un  — >•  v  in  C(D) .  In 
particular,  unj(x)  — >  v(x),  so  u(x)  =  v(x).  Since  un  >  u  in  D  for  each  n  6  N,  it 
follows  that  v  >  u  in  D.  □ 

Lemma  2.8  J"  =  Winf 

PROOF:  Let  u  =  |  hmu„  (E  J,  un  £  ftinf.  Fix  i£D.  By  Lemma  2.7,  for  each 
n  G  N  there  is  vn  e  H  such  that  i>n  >  un  in  D  and  vn(x)  =  «n(x).  Hence  {vn} 
is  bounded  at  x  and  therefore  relatively  compact.  Thus,  there  exists  a  subsequence 
{vnj}  and  vx  eH  such  that  vn}  -►  v*  in  C(D) .  Furthermore, 

vx(x)  =  \imvnj(x)  =  limunj(x)  =  u(x). 
Since  vnj  >  uUj  in  C(D) ,  it  follows  that  v*  >  u  in  Z). 

Let  Ru  =  ini{v  €  ft  :  v  >  u}.  Then  Ru  <E  ftinf  and  #u  >  u.  But,  for  a:  G  79,  vr  >  u 
and  v*(x)  =  u(x).  Hence,  Ru(x)  <  u(x).  This  proves  Ru  =  u,  i.e.,  u  6  Hini.  a 
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By  putting  together  previous  lemmas,  we  obtain  the  following. 

THEOREM  2.1  Let  H  be  a  closed  convex  cone  of  strictly  -positive  (unless  zero)  con- 
tinuous functions.  Assume  that  7i  has  a  compact  basis.   Then  the  family 

7iinI  =  {u  :  u  =  inf  uQ,  ua  £  H} 

is  a  closed  convex  cone,  stable  for  arbitrary  infima  and  closed  in  the  topology  of 
pointwise  convergence. 

Remarks:  (i)  Let  Tip  =  {u  £  "H  :  u(y)  =  1},  y  £  D.  Then  7{y  is  a  compact  basis  for 
H.  Therefore,  by  (2.2), 

u(  X ) 
k(x,y)  =    inf   — —  =  inf{u(:r)  :  u  G  H,  u(y)  =  1}. 

Let  us  denote  the  function  above  by  sy(x).  This  function  evaluated  at  x  gives  the 
greatest  lower  bound  on  numbers  u(x)  for  all  u  6  7Y,  given  that  u(y)  =  1.  Moreover, 
this  bound  is  attained  in  the  following  sense:  For  given  x  6  D,  there  exists  v  £  H 
(depending  on  x  and  y)  such  that  v  >  sy  and  v(x)  =  sy(x)  (see  Lemma  2.7). 
(ii)  W0*  is  a  closed  convex  cone,  and  (Hint)y  is  a  compact  basis.  One  can  consider 
the  function 

k(x,y)=     inf-—. 

u6?fnf  u(y) 

But  this  function  is  equal  to  k(x,y),  so  nothing  new  is  achieved.  Note  that  H"31  also 
satisfies  Harnack  principle. 

In  the  sequel  we  consider  some  other  properties  of  the  cone  H"**.    First  we 
establish  a  large  class  of  functions  which  operate  on  Hild . 

PROPOSITION  2.2  Assume  that  1  £  Hmi .  If<f>  :  (0,  oo)  -►  R  is  a  positive,  increasing, 
concave  function,  then  <f>  o  u  £  Hild  for  each  u  £  HinI . 


IS 


PROOF:  <j>  is  an  infimum  of  linear  functions  l(t)  =  at  +  b,    a  >  0,  b  >  0.  If  it  G  WiBf, 
then  /  o  u  =  au  +  b  is  again  in  Hlxd.  Let  ^  =  inf  (,-.  Then 

^oti  =  (inf  /,)  o  u  =  inf  (L  o  u)  €  H". 

D 

Another  operation  which  leaves  7{mi  invariant  is  integration.  Let  E  be  a  com- 
pact topological  space  with  the  Borel  <7-algebra  B.  Assume  that  u:DxE->  R+ 
has  the  following  two  properties: 

(i)  x  i-»  u(x,  <t)  G  "Winf  for  each  o  £  E, 

(ii)  a  i-*  u(x,er)  is  continuous  for  each  x  £  D. 

Let  ^  be  a  positive  Radon  measure  in  (E,  5).  By  Theorem  12.11  in  [6],  the  cone 
generated  by  point  mass  measures  {e^},  a  6  E,  is  dense  in  the  cone  of  all  positive 
Radon  measures  A^+(E)  (in  vague  topology).  Hence,  there  is  a  sequence 
{/in}  C  M+(YJ)  of  measures  of  the  type 

k 

l*n  =  '52aie<7i,    a,  >  0, 

which  converges  to  \i.  Thus,  for  every  /  £  C(E)  , 

/  /  dfin  ->        f  dp. 
In  particular,  for  each  x  6  D, 

/  u(x,a)nn(dcr)    -*    /  u(x,a)fi(dcr). 
But, 

f  k 

/  u(x,(r)fin(d(T)  =  J2aMx,°)  e  ftinf, 
"/E  .=1 
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and  so  /E  u(x,  a)fi(da)  G  Winf ,  because  7iinf  is  closed  for  pointwise  convergence. 

Let  E  be  a  locally  compact  space  with  countable  basis,  and  let  {£„}  be  an 
increasing  sequence  of  compact  subsets  of  £  such  that  E  =  UnSn.  If  /z  is  a  positive 
Radon  measure  on  the  Borel  cr-algebra,  then  for  fin  —  /j,\zn  we  have  that 

/     u(x,a)nn(d(T)  =   /  l^nu(x,a)fi(da)  G  Hliai. 

The  Lebesgue  monotone  convergence  theorem  and  pointwise  closedness  of  7Yinf  imply 
that  f^u(x,a)fi(d(r)  G  Hmf  (unless  identically  +00). 

Let  us  assume  now  that  the  cone  Ti  contains  positive  constants  and  separates 
points.  The  family  T>  =  Winf  —  7{ini  is  a  vector  lattice  of  continuous  functions  in  D, 
containing  all  constants  and  separating  points.  By  the  Stone- Weierstrass  theorem, 
T>  is  dense  in  C(D) .  Moreover,  if  /  €  C(D)  is  bounded  by  M,  and  un  — ►  /  where 
un  G  V,  then  vn  =  (un  A  M)  V  (-M )  G  V  and  vn  -»■  /  in  C(£>)  boundedly. 

Therefore,  if  ^x  is  a  positive  Radon  measure  on  the  Borel  <7-algebra  on  D  and 
/  G  C(-D)  with  compact  support,  then  the  Lebesgue  convergence  theorem  implies 
that  vn  — >  /  in  CP(D,  fi),  1  <  p  <  00,  where  {vn}  is  a  bounded  sequence  of  functions 
in  T>  converging  to  /  in  C(D) .  Since  continuous  functions  with  compact  support  are 
dense  in  £P(Z),//),  it  follows  that  V  is  dense  in  Cp(D,fi). 

For  the  next  consideration  we  assume  that  the  cone  H  is  linearly  separating: 
For  x,y  G  D  and  A  >  0,  there  exists  u  G  H  such  that  u(x)  ^  Xu(y). 

Let  Cq(D)  denote  the  set  of  continuous  functions  that  vanish  at  infinity.  We 
shall  assume  that  the  family  Hjf  =  HiT*  D  C0(D)  is  nonempty.  Let  u  G  H^  and> 
without  loss  of  generality,  we  assume  that  u(x0)  =  1  where  x0  is  some  fixed  point  in 
D.  For  given  e,  let  K  be  a  compact  set  such  that  u(x)  <  t  in  D  \  K.  Then 

c  >  u(x)  >  k(x,x0)u(x0)  =  sxo(x) 
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for  x  G  D  \  K,  and  so  sxo  G  Tl.^  ■  From  inequality  (2.4)  it  follows  that  sy(x)  < 
sXo(x)/sXo(y)  and,  therefore,  sy  €  H1^   for  every  y  G  D. 

Note  that  WJJ^  is  also  a  convex  cone  stable  for  arbitrary  infima.  If  u  G  "Hinf 
and  t>  is  an  element  of  H^  ,  then  un  =  (nu)  A  u  G  WJJ^  ,  for  each  n  £  N,  and  un  |  u. 
Hence,  every  function  in  7{inI  is  an  increasing  limit  of  functions  in  Ti1^1  .  In  particular, 
Ho*   is  also  linearly  separating. 

Following  [2]  we  say  that  a  convex  cone  V  C  C(D)  is  a  function  cone  if  V 
satisfies  the  following  three  conditions: 

(Fi)  There  exists  v  G  V  such  that  v  >  0; 

(F2)  V+  is  linearly  separating; 

(F3)  "P  is  adapted,  i.e.,  for  each  v  G  "P  there  exist  «GP  such  that  for  each  e  >  0, 

there  exists  a  compact  set  K  in  Z)  with  the  property  that  v(x)  <  eu(x)  for  every 

xeD\K. 

Lemma  2.9  W^1    is  a  function  cone. 

PROOF:  Both  (Fi)  and  (F2)  follow  from  assumptions.  For  v  G  Kg*  ,  let  u  =  y/v 
and,  for  e  >  0,  let  /(  =  {i  6  D  :  y/v(x)  ^  cl-  From  Proposition  2.2,  u  G  Winf, 
and  so  u  G  Wj)nf  .  Since  y/v  G  C0(I>)  ,  K  is  compact.  Further,  for  x  G  D  \  if, 
v(x)  <  eyv(x)  =  eu(:r)  which  proves  (F3).  n 

It  is  not  difficult  to  see  that  7iinf  need  not  be  a  function  cone.  One  explicit 
example  will  be  given  in  Chapter  4. 

As  it  should  be,  reasonable  functionals  on  Ti"^  are  given  by  measures.   Let  $ 
be  an  additive,  positively  homogeneous,  increasing  functional  on  Wini.  Then  $  n,ini 
has  the  same  properties.   By  one  of  the  numerous  Choquet's  results,  there  exists  a 
unique  measure  \i  on  the  Borel  cr-algebra  of  D,  such  that  $(u)  -  fDv  dfi  for  every 
v  G  Hf   (e.g.,  [2,  p.17,  Prop.1.4]). 
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If  u  G  "Hinf,  let  {un}  be  an  increasing  sequence  of  functions  in  7^q  such  that 
u  =t  limn  un.  Then 

$(u)  =  lim$(un)  =  lim  /   un  da  =   /   u  o?u. 

This  proves  the  following. 

PROPOSITION  2.3  Let  $  :  7iint  — ►  R  be  an  additive,  positively  homogeneous,  increas- 
ing functional.  Then  there  exists  a  unique  measure  n  on  D  such  that  $(u)  =  / u  du 
for  every  u  6  'Had . 

A  function  u  €  'Hmi  is  said  to  be  extremal  if  u  =  v  +  w  for  u,  w  €  7ilnf  implies 
that  u  and  u;  are  proportional  to  it,  i.e.,  there  is  A  >  0  such  that  v  =  Au,  w  =  (1  —  A)u. 
We  show  that  for  each  y  (z  D  the  function  Sj,  =  k(-,y)  is  extremal  in  7iinf. 

If  sy  =  v  +  w  where  v,w  €  'Hinf,  then  1  =  Sj,(?/)  =  v(y)  +  w(y).  Further,  for 
x  G  Z?,  u(a;)  >  k(x,y)v(y)  and  it>(ar)  >  k(x,y)w(y).  Therefore, 

sy(a:)  =  v(x)  +  w(x)  >  sy(x)v(y)  +  sy(x)w(y)  =  sy(x). 

Hence,  v(x)  —  v(y)sy(x)  and  w(x)  =  w(y)sy(x).  Since  this  holds  for  each  x  G  D,  we 
have  u  =  v(y)sy  and  t«  =  w(y)sy. 

We  will  show  in  Chapter  4  that,  in  general,  there  are  more  extremal  functions 
than  just  sy,  y  G  D. 

2.2     Kernel  Function 

Let  D  be  a  topological  space  having  the  same  properties  as  in  the  previous 
section.  In  addition,  we  assume  that  D  is  contained  in  a  compact  metrizable  space 
denoted  by  D,  such  that  D  is  the  interior  of  D  and  the  metric  of  ~D  restricted  to  D 
is  d.  Let  dD  denote  D  \  D;  we  call  dD  the  boundary  of  D. 
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Let  H  be  a  closed  convex  cone  in  C+(D)  containing  the  function  identically 
equal  to  zero.  We  also  assume  that  there  is  u0  6  H  satisfying  m  <  u0  <  M  for  some 
strictly  positive  constants  m  and  M.  Let  x0  be  an  arbitrary,  but  fixed  point  in  D. 

In  this  section  we  shall  assume  the  existence  of  a  function  on  D  x  dD,  which 
we  call  the  kernel  function.  In  the  next  chapter,  it  will  be  shown  that  such  function 
exists  for  a  broad  class  of  convex  cones. 

The  basic  hypothesis  is 

(Hi)  There  exists  a  function  K  :  D  x  dD  — >  R  such  that 

(i)  for  each  z  (E  dD,  x  t->  K(x,z)  belongs  to  Hxo., 
(ii)  for  each  x  E  D,  z  h->  K(x,z)  is  continuous  on  dD. 

Before  we  give  the  second  hypothesis,  let  us  note  the  following  simple  fact:  If 
fi  is  a  finite  measure  on  the  Borel  subsets  of  dD,  then  the  function  u  defined  by 


u(x)  =  /     K(x,z)fi(dz) 

JdD 


is  in  H.  This  follows  from  one  of  the  results  of  Section  2.1  (see  p. 19).  We  shall  assume 
that  all  functions  in  7i  arise  in  this  way. 

(#2)  For  each  u  £?{,  there  exists  a  unique  Borel  measure  fi  on  dD  such  that 


u(x)  =  /     K(x,z)n(dz). 

Jon 


IdD 

Note:  Since  K(x0,z)  =  1  it  follows  that  u(x0)  =  fi(dD),  so  /x  is  a  finite  measure. 

Let  us  recall  the  notation  from  the  previous  section:    For  a  6   D,  sa(x)  = 
inf{u(x)  :  u  G  H,u(a)  =  1}.  Let 

Ha,x  =  {ueH:  u(a)  =  1,  u(x)  =  sa(x)},    a,x  <E  D.  (2.12) 
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From  the  Remark  following  Theorem  2.1,  it  follows  that  7ia<x  is  nonempty.  It  is  also 
convex,  compact  and  closed  for  pointwise  convergence.  Since  every  function  in  Ti  is 
representable  by  the  kernel  function,  it  is  reasonable  to  expect  that  7iaiX  contains 
functions  of  the  form  x  h-»  K(x,z),  z  6  dD  (properly  normalized).  We  show  that 
this  is  true. 

Let  u  6  Tia,x  aQd  H  the  measure  representing  u.  Then 

sa(x)  =  u(x)=  f    K(x,z)(i(dz).  (2.13) 

JdD 

Assume  that  A'  and  A"  are  disjoint  Borel  subsets  of  dD,  such  that  A'  U  A"  =  dD 
and  fi(A')  >  0,  /x(A")  >  0.  Let  a'  =  fA,  K(a,z)y,(dz)  and  a"  =  fA„  K(a,z)p(dz). 
Then  both  a'  and  a"  are  strictly  positive,  so  one  can  define  measures  v'  and  v"  on 
dDby 

v  =  —fJ,\A'    and   v    =  — j,[i\tj>- 


Let 


u'=  I    K(-,zy{dz)   and   u"  =  I    K(;z)v"(dz). 

JdD  JdD 


IdD 

Then  u'  and  u"  are  in  H  and  a  simple  computation  shows  that  u'{a)  =  u"{a)  =  1. 
Hence,  u'(x)  >  sa(x)  and  u"(x)  >  sa(x).  Therefore,  since  a'  +  a"  =  1, 

u(x)  =  a'u'(x)  +  a"u"{x)  >  a'sa(x)  +  a"sa(x)  =  sa(x)  =  u(x). 

This  implies  that  u'(x)  =  u"(x)  =  sa(x).  Hence,  u'  and  u"  are  in  H^x  and  represent- 
ing measures  have  smaller  support. 

Proposition  2.4  For  a  G  D  and  x  e  D,  there  exists  z  =  z(a,  x)  €  dD  such  that 

KM    e   H 
K(a,z) 
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PROOF:  Let  u  eHax  and  u  =  /     K(-,  ()fi(d().  If  fi  is  a  multiple  of  a  point  mass  at 

JdD 

z,  there  is  nothing  to  prove.  If  //  charges  some  point  z  €  dD,  we  take  A'  =  {z}  in  the 
construction  preceding  the  statement.  Then  the  function  u'  from  above  is  precisely 
K(-,z)/K(a,z). 

So  we  assume  that  \i  does  not  charge  points.  Let  z  G  Supp//.  Then  //  charges 
every  neighborhood  of  z  in  dD.  Let  {An}  be  a  decreasing  sequence  of  neighborhoods 
of  z  in  dD  which  shrink  to  z,  i.e.,  fl^LjAn  =  {z}.  There  is  n  €  N  such  that 
fi(dD  \  An)  >  0,  since  otherwise  \i  would  be  concentrated  at  z.  We  assume,  without 
loss  of  generality,  that  this  is  true  for  every  n  (E  N. 

Let 

an  =  /     K(a,()fi(d(),     vn= — //|An     and     un  =  \     K(-,Qvn(d(). 
•/a„  an  JdD 

Then  «n  6  7iax  for  each  n  €  N. 

We  claim  that  the  sequence  of  measures  {vn}  is  bounded.  First  note  that 

„„m .  e&J .    //A;>     .  (2.i4) 

Since  £  *-*  A^a,  ()  is  continuous  on  dD,  there  is  n0  G  N  such  that  for  every  £  G  A,^ 

iff(a>*)<ff(a,C)<2Jir(a,*). 
Integrate  the  above  inequalities  with  respect  to  //  over  An,  n  >  n0,  to  get 

±ff(a,z),,(An)  <  /    K(a,()ti(dQ  <  2K(a,zMAn). 
After  dividing  ^(An)  by  each  of  the  terms  above,  we  obtain 

1  <  MAn)  2 

2A'(a,z)  -  /An  A'(«,C^K)  "  K(a,z)  ' 
But  the  middle  term  is,  by  (2.14),  equal  to  vn(dD),  which  proves  the  claim. 
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Boundedness  of  the  sequence  {vn}  implies  the  existence  of  a  weakly  convergent 
subsequence.  We  may  assume,  without  loss  of  generality,  that  {un}  weakly  converges 
to  a  positive  Borel  measure  v  on  3D.  Let 

v=l    K(;Qv{d()  =  \imt    ff(.,CK(«*C)  =  li?i«n. 

JdD  n     JdD  n 

Since  'Ha,!  is  closed  for  pointwise  convergence,  v  £  H1:I. 

We  claim  that  v  is  concentrated  on  {z}. 

For  n  >  2,  let  /  be  a  continuous  function  on  dD  such  that  /  =  1  on  dD  \  An_i, 
/  =  0  on  ClAn_i  C  An  and  0  <  /  <  1.  Then,  for  every  m  >  n,  /  =  0  on  Am,  and, 
therefore,  jdD  fiQv^dQ  =  0.  But 

JdD  JdD 

Hence, 

0  =  L /(0"W)  >  /m A     HCHdO  -  K«U  \  An_a). 

•'3D  ^9£>\An_i 

Thus, 

K5Z)  \  {z})  =  i/(U~  !(^  \  A„))  =  0, 

so  v  sits  at  {z}. 

Therefore,  v  =  ctz  for  some  positive  constant  c.  Further 

1  =  ttn(fl)  =  [    K{a,Qvn(dC)  -  /    K(a,C)HdQ  =  cK{a,z). 

JdD  JdD 

K(-  z\ 

Hence,  v  =   ..;''    {  £  Wa>SB.  a 

K{a,z) 

It  is  evident  that  the  kernel  function  K  on  D  x  dD  has  the  central  role  in  all 
considerations  about  the  cone  H.  On  the  other  hand,  we  saw  in  Section  2.1  that  the 
function  k(x,y)  =  ini{u(x)/u(y)  :  U  £  H)  (denoted  by  sy(x)  here)  had  distinguished 
place  among  the  functions  in  HinI.  It  would  be  nice  if  the  kernel  function  K  could 
be  recovered  from  the  function  k.  This  is  true  with  the  additional  hypothesis  on  K : 
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(H3)  For  all  zx,z2  E  dD  such  that  zx  ^  z2, 

,        xlijP  v^(X>Z)    =0' 

where  (x,  z)  E  D  x  dD  and  (x,z)  — »  (21,22)  in  D  x  D. 

This  hypothesis  is  in  some  sense  stronger  than  the  previous  two  and  will  restrict 
the  examples  in  the  next  chapter  to  nice  domains.  Still,  as  we  shall  see,  there  are 
many  interesting  examples  satisfying  (H3). 

Let  us  fix  x  E  D  {xq  is  still  fixed).  For  each  a  E  D,  there  exists  z  —  z(a,  x)  E  dD 
such  that  K(-,z)/K(a,z)  E  7ia,x-  Denote  z(a,x)  by  za. 

Recall  that  we  have  assumed  the  existence  of  the  function  uQ  E  ti  satisfying 
m  5;  Uq  <  M.  For  any  a  in  D,  u0(x)/u0(a)  <  M/m.  Therefore,  sa(x)  <  M/m. 

LEMMA  2.10  Let  a  —>  z,  z  e  dD.  Then  za  -»  z. 

PROOF:  Let  {an}  be  a  sequence  in  D  converging  to  z  and  let  us  denote  the  corre- 
sponding points  on  the  boundary  by  {zn}.  Since  dD  is  compact,  we  may  assume  that 
{zn}  converges  to  some  point  z0  E  dD.  If  z0  ^  z,  then  by  (i/3),  limn  K(an,zn)  =  0. 
By  the  continuity  of  K,  limn  K(x,  zn)  =  K(x,zQ)  <  00.  Therefore,  the  sequence 
{K(x,  zn)/K(an,zn)}  is  unbounded.  On  the  other  hand, 

K(x,  zn)  M 

-777 r  =  san{x)  <  — . 

K{an,zn)  m 

Contradiction!  Hence  z0  —  z.  D 
Theorem  2.2  For  every  z  e  dD, 


lim      /  >  =  K(-,z). 
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PROOF:  Recall  that  x0  £  D  such  that  K(x0,  •)  =  1  on  3D.  Fix  x  E  D.  Let  va  and  za 
be  points  on  dD  such  that  K(-,va)/K(a,va)  €  %<,,*<,  and  K(-,za)/K(a,za)  e  "Ha,*, 
a  E  D.  Thus, 


=  sa(x)   and     T^,_      N    =  Mxo)- 


By  definition  of  sa, 


K(a,  za) 


K(x,va) 


>  sa(x)   and 


K(a,va) 
From  (2.15)  and  (2.16)  it  follows 

K(a,za)       K(x,za) 


K(a,va) 

K(x0,za) 
K(a,  za) 


(2.15) 


>  sa(x0). 


(2.16) 


> 


and    K&l£  <  l 


(2.17) 


K(a,va)  ~  K(x,va)  K(a,va) 

As  a  — »•  z,  Lemma  2.10  gives  that  za  — ►  z  and  va  — >  z.    The  first  inequality 
above  and  continuity  of  K  give 

liminf  |1^4  >  liminf  §^4  =  |H  =  1. 
»-**    /qa,ua)         «-*    #(z,u0)       /sT(x,z) 

From  the  second  inequality  in  (2.17)  it  follows  that 


Hence, 


By  (2.15), 


7i(a,za) 

hrnsup  — <  1. 

a—-z     K(a,va) 


BmStil-l. 


•-»  #(a,  ua) 


(2.18) 


sa(x) 
Sa(aro) 


Therefore,  by  (2.18)  and  continuity  of  K, 


K{x,  za) 

K(a,za) 

K(x0,va) 
K(a,va)  . 

A"(ar,  Z„) 


K(a,za)' 


lim  — r-^-  =  hm  K(x,  za)  Ty.)_  _  [  =  K(x,  z) 


a^z  sa(x0)       ■»-**' 


K(a,za) 
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CHAPTER  3 
EXAMPLES  OF  CONES  WITH  COMPACT  BASIS 


3.1     Positive  Solutions  of  Elliptic  Differential  Equations 
Let  D  be  a  bounded  domain  in  Rn,  and  let 

«'J=i 
be  an  elliptic  operator  in  divergence  form  with  bounded  coefficients  <z,j  satisfying 
a'j  —  aji-  We  shall  assume  that  A  is  uniformly  elliptic  in  D,  i.e.,  there  exist  constants 
A  and  A,  0  <  A  <  A,  such  that 

AKI3<£M*)tti<A|£|2 

for  all  £  G  Rn  and  x  €  D. 

A  weak  solution  of  the  equation  Au  =  0  in  D  is  a  function  u  in  the  Sobolev 
space  Htfc(D)  =  {ue  qoc(D)  :  Vu  €  Lfoc(D)}  satisfying 

n 

^    /   aij(x)Diu(x)Dj<f>(x)  dx   =  0  (3.1) 

for  every  (f>  G  Cf{D). 

It  is  well-known  that  every  solution  of  Au  =  0  is  locally  Holder  continuous  in 
D:  For  each  compact  subset  K  of  D 

\u(x)  —  u(y)\ 

sup  J-V-i -^±  <  oo,    where  0  <  a  <  1. 

x,yeK      \x  -  y\a 

This  is  a  famous  result  discovered  by  DeGiorgi  and  Nash.  See  [13,  p.  190]  for  a  proof. 
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In  particular,  local  Holder  continuity  of  solutions  has  as  a  consequence  the  fact 
that  solutions  are  locally  uniformly  equicontinuous  (and,  of  course,  continuous  on  D). 

The  next  result  we  need  is  Harnack  inequality  proved  by  Moser.  A  proof  can 
be  found  in  [13].  A  more  convenient  form  is  the  one  given  in  [4].  Let  K  be  a  compact 
subset  of  D.  Then  there  exists  a  constant  c  =  c(A,  A,K,  D,n)  such  that  for  every 
nonnegative  solution  of  Au  =  0, 

maxu  <  crmnu. 
K        ~       K 

By  linearity  of  A,  nonnegative  solutions  of  Au  =  0  form  a  convex  cone.  This 
cone  is  the  object  of  our  interest.  Formally,  let 

H=  {u:  Au  =  0,u  >  0}.  (3.2) 

All  functions  in  7i  are  continuous.  If  u(x)  =  0  for  some  x  in  D,  u  €  7i:  then  Harnack 
inequality  immediately  implies  that  it  is  identically  zero.  Thus,  H  consists  of  strictly 
positive  functions  (unless  zero). 

Let  us  show  that  Tt  is  closed  in  C(D) .  The  following  estimate  is  needed. 

LEMMA  3.1  Let  y  G  D  and  0  <  r  <  R  such  that  the  ball  B(y,  R)  is  contained  in  D. 
For  any  solution  u  of  Au  =  0  the  inequality 

I        |Vu(x)|2  dx<C  I         u2(x)  dx  (3.3) 

JB(v,t)  JB(y<R)  v       ' 

holds,  where  the  constant  C  depends  on  A,  A,r  and  R. 

This  lemma  is  standard  and  is  stated  in  a  similar  form  in  [18].  For  a  proof, 
first  note  that  (3.1)  is  also  valid  for  <j>  G  H~l'2(D).  Let  V  be  C~  function  identically 
1  in  B(y,r),  zero  outside  B(y,R)  and  with  bounded  gradient.  Let  <f>  =  ux/>2;  then 
(f>  G  Hq   .  Equation  (3.1)  can  be  written  as 
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—  /         ib2y^  ciijDiuDjU  dx  —  2  /         tz0  y  D{uDjib  dx. 

JB(y,R)        *rf  JB(y,R)         ** 

Repeated  use  of  ellipticity  and  Schwartz  inequality  proves  the  lemma. 

Now  let  {uk}  be  a  sequence  in  7i  and  assume  that  uk  —*  u  in  C(D) .  First  we 
show  that  u  is  in  H* '  (D).  Let  y€.D,0<r<R  with  B(y,  R)  C  -D.  Since  uk  —  wm 
is  a  solution  of  Au  =  0,  the  lemma  above  yields 

/         |V(ufc-um)|2  dx  <  C  /         (uk-um)2  dx. 

JB(y,r)  JB(y,R)y 

Since  Ujt  — ►  «  uniformly  in  B(y,  R),  the  integral  on  the  right  hand  side  is  arbitrarily 
small  for  large  k  and  m.  Hence,  the  left-hand  side  is  also  small.  Therefore,  for  t  >  0, 
there  is  k0  G  N  such  that  for  k,m  >  k0 


I        (\uk  -  um\2  +  |V(ufc  -  um)\2)dx  =||  uk  -  um  \\2Hi,2(B(ytr)) 


<  e. 


Here  H1,2(B(y,r))  denotes  Sobolev  space  on  B(y,r).  By  completeness,  there  exists 
v  £  H1'2(B(y,r))  such  that  uk  ->  v  in  Hh2(B(y,r)).  In  particular,  uk  — >  u  in 
L2(B{y,r)).  But  ufc  — ►  u  uniformly  on  B(y,r),  so  u  =  v  a.e.  in  B(y,r).  Therefore, 
«  G  H1'2(B(y,r)).  Since  y  and  r  were  arbitrary,  u  E  fli1"3  (£>). 

Furthermore,  uk  — ►  u  in  H1'2(B(y,r))1  which  implies  that  u*.  -►  ix  in  Hia(Dx) 
for  every  relatively  compact  open  subset  Di  of  D. 

If  <£  G  C~(Z>)  with  Supp(<£)  C  Z?j,  then  (3.1)  holds  for  every  uk.  Note  that 
maxj  ||  Dj(j)  ||oo<  oo  and  max.j  ||  a^  ||o0<  oo.  Let  c  be  a  common  bound  for  both 
maxima.  Then 

JD  \J2aiAx)(D'uk  ~  Diu)Dj<j)\  dx 
<  JDi  [E  II  «0  He.  \DiUk  -  DiU\  ||  Drf  HooJ  dx  <  c2  f^  (j2  \DiUk  -  DM  dx 
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=  nc2J    (Y,\Diuk-Diu\\dx<nc2Yjrn{DxY'2U    |AK  -  u)\2dx\ 

=  nc2  miD^'2  £  Qk  \Dt(uk  -  u)\2dx)  "  , 
which  tends  to  zero  since  uk  —»  u  in  H1,2(D\).  Thus 

52  /   a,ij{x)DiuDj4>  dx  =  y^      aij{x)(DiU  —  DiUk)Dj<f>  dx 

is  arbitrarily  small  for  large  A;,  hence  zero.    Since  <f>  was  an  arbitrary  function  in 
C£°(.D),  it  follows  that  Au  =  0.  Therefore,  u  G  7Y;  so  ?i  is  closed. 
Let  #0  be  an  arbitrary,  but  fixed,  point  in  D  and  let 


HXQ  =  {u:  u  e  7i,u(x0)  =  1}. 


(3.4) 


Then  7iXQ  is  closed  in  C(D) ,  locally  uniformly  equicontinuous  and  bounded  (by  Har- 
nack  inequality).  Further,  TiXo  is  obviously  a  basis  for  H.  Hence,  H  has  a  compact 
basis.  Note  that  Ti  contains  positive  constants.  This  shows  that  7{  is  a  cone  of  the 
type  studied  in  Chapter  2,  so  all  results  from  Section  2.1  are  applicable.  In  partic- 
ular, the  cone  Ti*"*  of  all  possible  infima  of  positive  solutions  of  Au  —  0  consists  of 
continuous  functions  and  is  closed  for  pointwise  convergence.  Of  course,  the  case  of 
the  Laplacian  A  =  £d2/dx-2  is  also  included. 

In  order  to  apply  results  of  Section  2.2,  some  restrictions  on  the  domain  D  are 
needed.  We  shall  assume  that  D  is  a  bounded  Lipschitz  domain  in  Rn. 

A  bounded  domain  D  C  Rn  is  called  a  Lipschitz  domain  if  there  are  positive 
numbers  r0  and  m  such  that  for  every  point  z  £  3D 

B{z,r0)  DD  =  B(z,r0)  n  {(x,t)  :  x  €  Rn~\t  >  <f>(x)},    and 

B(z,  r0)  D  3D  =  B{z,  r0)  (1  {(x,  <f>(x))  :  x  €  R""1}, 
where  (f>  is  a  function  satisfying  ||  V(j>  ||Loort>"-^<  m  <  oo. 
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We  denote  B(z,r)  D  dD  by  A(z,r).  For  s  >  1,  let 

^(z^sr)  =  {(y,t)  :  |y  -  y0|  <  r,\t  -  <f>(y0)\  <  sr},    z  =  (y0,  <f>(y0)), 

be  a  cylinder,  and  let  vr  =  vr(z)  =  (y0,  <f>(yo)  +  r)  be  the  "center"  of  that  cylinder. 

Let  us  recall  the  definition  of  the  elliptic-harmonic  measure.  The  Dirichlet 
problem  for  A  in  D  is  solvable  if  for  every  continuous  function  g  on  dD,  there  is  a 
unique  continuous  function  u  on  D  such  that  Au  —  0  in  D  and  u\g£>  =  g.  A  result 
from  [18]  states  that  the  Dirichlet  problem  for  A  in  D  is  solvable  if  and  only  if  it 
is  solvable  for  the  Laplacian  A.  Since  D  is  Lipschitz,  it  follows  that  the  Dirichlet 
problem  for  A  in  D  is  solvable. 

For  x  €  D,  the  mapping  g  \- >  u(x)  is  a  positive  linear  functional  on  C(dD). 
Hence,  there  is  a  positive  measure  lox  on  dD  such  that 


t(x)=  [    g(()u*(dO. 

J  3D 


It  is  easy  to  see  that  ux  is  a  probability  measure.   We  call  u>x  the  elliptic-harmonic 
measure  for  A  in  D  at  x. 

Fix  z  6  dD.   A  function  K(-,z)  defined  on  D  is  called  a  kernel  function  at  z 
(for  the  operator  A  and  domain  D)  if 

(i)  K(-,  z)  is  a  solution  of  Au  =  0  in  D, 

(ii)  K(;  z)  e  <:(£  \  {z})  and  limx^c  A'(x,  z)  =  0  for  (  e  dD,  C,  ±  z,  and 

(iii)  K(x,  z)  >  0  for  each  x  £  D  and  A'(x0,  z)  =  1,  where  x0  is  the  fixed  point  in  Z>. 

The  existence  of  the  kernel  function  for  A  is  established  in  [4,  Theorem  3.1]  with  the 
additional  property  that  the  function  z  i-+  K(x,z)  is  continuous  on  the  boundary 
dD  ([4,  Corollary  3.2]).  Therefore,  K,  considered  as  a  function  of  two  variables  on 
D  x  dD,  satisfies  the  hypothesis  (Hi)  from  Section  2.2. 
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The  representation  theorem  was  also  obtained  in  [4,  Theorem  4.1]  (see  also  [9]): 
For  each  u  G  7~t,  there  is  a  unique  Borel  measure  /j.  on  dD  such  that 


u(x)  =  /     K(x,z)[i(dz). 

JdD 


Thus,  our  hypothesis  (H2)  also  holds. 

To  establish  (H3)  we  need  to  show  that  whenever  (x,  z)  G  D  x  dD  converges  to 
(zi,z2),  Z\,z2  G  dD,  Z\  ^  z2,  the  kernel  function  K(x,  z)  tends  to  zero.  Once  again 
we  reach  for  a  result  from  [4].  This  time  it  is  their  Lemma  2.5 

LEMMA  3.2  Suppose  r  <  ro/2  and  u  is  a  positive  solution  of  Au  =  0  with 


(i)  u  G  C(D  \  tp(z2,  r,  sr)),  and 


(ii)  u  =  0  ondD\A(z2,r). 

Then,  for  x  G  D  \  ip(z2,2r,2sr), 

c-1u(ur)u;*(A(z2,r))  <  u(x)  <  cu(vr)ujx(A(z2,r)) 

(where  vT  =  vr(z2))  with  c  depending  only  on  A,  A  and  M. 

We  use  this  lemma  for  functions  x  \~*  K(x,z).  For  zx  ^  z2,  there  exists  r  >  0 
such  that  zx  0  ij>(z2,3r,3sr).  If  z  G  ip(z2,r,sr),  then  K(-,z)  satisfies  assumptions  (i) 
and  (ii).  Hence,  for  x  G  ip(z2,2r,  2sr), 

K(x,z)  <  cK(vr,z)u*(A(z2,r)). 

As  x  ->  zu  ujx(A(z2,r))  ->  0,  and  as  z  -►  z2,  K(vr,z)  ->  K(vr,z2)  <  oo. 
Hence,  K(x,z)  -+  0  as  (x,z)  ->  (z1,z2).  This  proves  (//3). 
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3.2     Positive  Solutions  of  Schrodinger  Equation 

In  this  section  we  consider  positive  solutions  of  the  Schrodinger  equation.  The 
notation  remains  as  in  Section  3.1:  D  is  a  bounded  domain  in  Rn,  but  here  we  require 
that  n  >  3. 

Let  L  =  —A  +  q  be  the  Schrodinger  operator  on  D.  A  is  a  uniformly  elliptic 
operator  in  divergence  form  as  in  3.1  with  bounds  A  and  A,  and  q  is  a  function  in  the 
Kato  class  Kn(D),  i.e., 

lim  sup  /  ,    l<?l(^  Jy  =  0.  (3.5) 

A  weak  solution  of  Lu  =  0  is  a  function  u  in  the  Sobolev  space  Hi'(D) 
satisfying 

—  5")   /  a,ij(x)Diu(x)Dj<f>(x)  dx  —       q(x)u(x)<f>(x)  dx  (3.6) 

for  every  function  <f>  £  C™(D). 

Continuity  of  solutions  and  Harnack  inequality  for  nonnegative  solutions  was 
established  in  [5].  They  proved  the  following  two  theorems. 

CONTINUITY  THEOREM:  There  exists  a  nondecreasing  function  u(s)  depending  on 
A  and  A  such  that  u(s)  — *  0  as  s  — ►  0  and,  for  any  solution  of  Lu  =  0  in  D,  and  for 

any  ball  BT{y)  with  BzT{y)  C  D, 

(\^£  —  y\\ 
sup  \u\,    x  e  Br(y).  (3.7) 
r       J  B3r(y) 

HARNACK 's  THEOREM:  There  exist  positive  constants  ro  and  C  depending  on  X,  A 
and  n,  such  that  for  any  nonnegative  solution  u  of  Lu  —  0  in  D  and,  for  any  ball  BT 
with  0  <  r  <  r0  and  B2r  C  D,  we  have 

supu    <   C  inf  u. 

fir/2  B-72 
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From  Harnack's  Theorem  it  follows  that  for  every  compact  subset  K  of  D  there 
is  a  constant  C  =  C(\,A,n,K)  such  that  for  any  nonnegative  solution  u  of  Lu  =  0 
in  D  we  have 

supu   <   Cinfu.  (3.8) 

This  follows  from  the  familiar  chain  argument. 

Let  H.  denote  the  family  of  all  nonnegative  solutions  of  Lu  =  0  in  D.  Note  that 
Harnack  inequality  implies  that  every  nonnegative  solution  which  is  not  identically 
zero  is  necessarily  strictly  positive.  Let  xq  be  an  arbitrary,  fixed  point  in  D  and  let 

HXo  =  {u:u<EH,u(x0)  =  l}.  (3.9) 

Harnack  inequality  immediately  shows  that  HXo  is  bounded  in  C(D) .  Let  y  be  any 
point  in  Z),  and  choose  r  >  0  such  that  B3r(y)  C  D.  Then  H^  is  bounded  on  B3r(y). 
Now  (3.7)  implies  the  uniform  equicontinuity  of  HXo  at  y.  Since  this  holds  for  every 
y  G  D,  7ij.0  is  locally  uniformly  equicontinuous. 

It  remains  to  show  that  7i  is  closed  in  C(D) .  This  is  proved  in  the  exactly  same 
way  as  in  Section  3.1.  Lemma  3.1  is  valid  for  solutions  of  Lu  =  0,  and  the  proof  can 
be  found  in  [5]. 

Hence,  H  is  a  closed  convex  cone  with  compact  basis  HXo .  Thus,  we  can  form  the 
cone  Hint  consisting  of  all  possible  infima  of  positive  solutions  of  Lu  =  0.  Functions  in 
H™*  are  continuous,  and  Ww  is  closed  for  pointwise  convergence.  So  again,  everything 
worked  fine.  But,  contrary  to  the  case  of  the  operator  A,  this  time  H  may  consist  of 
the  zero  function  only. 

In  order  to  obtain  non-void  results,  both  the  domain  D  and  the  function  q  need 
to  be  specialized.  As  in  3.1,  we  suppose  that  D  is  a  bounded  Lipschitz  domain.  The 
function  q  must  have  finite  gauge.   To  see  what  that  means  we  follow  [9].   We  will 
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only  outline  the  set-up.  The  detailed  discussion  can  be  found  in  the  aforementioned 
paper. 

Let  dx  denote  the  Lebesgue  measure  on  D.  The  symmetric  form  £  on  L2(D,  dx) 
associated  with  A  is  defined  by 


£(u,v)  =  Y,   f  aij(x)Diu(x)Djv(x)  dx,  V[£]=C^(D). 

This  form  is  Markovian,  meaning  that  the  following  axiom  is  satisfied: 


(3.10) 


For  each  e  >  0  there  exists  a  real  function  <j)e(t),  t  G  R,  such  that 

<j>t{t)  =  t    Vie  [0,1],     -e<<f>e(t)<l  +  e    VI6R,  (3.11) 

and  0  <  <f>t(t')  -  <f>c(t)  <t'  -t   whenever  t  <  tf, 

u  e  V[£]  =>  <f>e(u)  G  T>\S\  and  £(&(«),  Mu))  <  £(u> u)  (3-12) 

(see  [12,  p.7]). 

The  form  (3.10)  has  the  local  property,  i.e.,  if  u,v  6  T>[£],  and  Supp(u)  and 
Supp(u)  are  disjoint  compact  sets,  then  £{u,  v)  —  0. 

The  form  (3.10)  is  not  closed,  but  it  is  closable  ([12,  p. 43]).  Its  smallest  exten- 
sion £  is  a  regular  Dirichlet  form  possessing  the  local  property  ([12,  pp. 41, 42]).  This 
means  that  £  is  a  closed,  symmetric  Markovian  form  possessing  a  core.  A  core  of  £ 
is  by  definition  a  subset  C  of  T>[£]  f\C0(D)  such  that  C  is  dense  in  T>\£]  with  respect 
to  the  norm  ||u||  =  £(u,  u)  +  (u,  u),  and  C  is  dense  in  Co(D)  in  the  uniform  norm.  By 
Theorem  6.2.2  in  [12],  there  exists  a  diffusion  (Xt,Px)  on  D,  with  no  killing  inside 
D,  whose  Dirichlet  form  is  the  one  given  by  £. 

We  assume  that  (Xt,Px)  is  defined  on  the  canonical  space  of  paths  (Q,Ft,0t). 
Here 

fi  =  {a;  :  t  i— ►  u(t)  is  continuous  on  (0,  ((w)),  w(t)  =  d,  t  >  C^)}) 
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where  ((uj)  is  the  lifetime  of  u>  and  d  denote  the  "cemetery".  The  filtration  (J7^  is 
natural  and  6t  is  the  shift. 

We  also  need  the  conditional  diffusion  starting  at  x  and  dying  at  y.  There  is  a 
measure  P*  on  f2  such  that  P*(X0+  =  x)  =  1  and  Py(XTD_  =  y)  =  1  where  Tjj  is  the 
killing  time  of  X  (i.e.,  the  exit  time  from  D).  For  ifD,  such  measure  is  constructed 
by  using  Doob's  A-transform,  while  for  x  6  dD,  time- reversal  is  also  used. 

For  the  function  q,  let 


e«v 


,(<)  =  exp{-  fq{Xa)ds}-  (3.13) 

Jo 

The  gauge  of  (A,  D,q)  is  defined  by 

F(x)  =  E*[eq(rD)],    xeD,  (3.14) 

and  the  conditional  gauge  by 

F(x,y)  =  Exy[eq(rD)],    x,y  eD.  (3.15) 

Here  tq  is  the  exit  time  from  D,  i.e.,  the  lifetime.  If  the  gauge  F  is  not  identically 
infinite,  then  it  is  bounded.  Moreover,  finiteness  of  the  gauge  is  equivalent  to  either 
of  the  following  two  conditions: 

(i)  There  is  a  solution  in  D  to  Lu  =  0  with  inf^  u  >  0. 

(ii)  For  each  /  G  C(dD)  with  /  >  0,  there  is  a  nonnegative  solution  to  the  Dirichlet 

problem 

Lu(x)  =  0,  x  e  D,      u(z)  =  /(z),  z  edD 

(for  the  proof  see  [9]). 

From  now  on  we  assume  that  F  ^  oo.  Then  the  cone  H  is  nontrivial.  Further- 
more, from  (ii)  it  follows  that  there  is  a  function  in  H  which  is  bounded  from  below 
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(by  a  strictly  positive  constant)  and  from  above.  Hence,  H.  satisfies  the  assumptions 
from  2.2.  It  remains  to  check  hypotheses  (i7i)-(^3). 

In  [9]  it  is  proved  that  the  conditional  gauge  F(-,  •)  is  continuous  on  D  x  D  and 
that  there  exist  constants  ca  and  c2  ,  0  <  C\  <  c2  <  oo,  such  that  cx  <  F(x,y)  <  c2, 
for  all  (x,y)  E^xB  (see  their  theorems  4.8  and  4.6).  By  Theorem  5.5  in  [9],  every 
function  u  (E  7i  has  the  unique  representation 

u(x)=  I    KL(x,z)fi{dx)  (3.16) 

JdD 

for  some  Borel  measure  /z  on  dD,  where 

Kl(x,z)  =  §^\k(x,z)  (3.17) 

(recall  that  xQ  is  fixed  in  D  and  K  is  the  kernel  function  for  A).  Continuity  of  KL  in 
both  variables  follows  from  (3.17)  and  the  continuity  of  K  and  F.  Finally, 

lim         KL(x,  z)<—  lim         K(x,  z)  =  0, 

(x,  z)—  (zi.za)  Ci       (x,z)^(zi,Zi) 

for  all  zi,z2  €  dD,  zi  /  ^2-  Hence,  all  three  hypotheses  are  satisfied. 

3.3     Harmonic  Functions  for  a  Diffusion 

In  this  section  we  switch  from  the  analytic  setting  to  a  probabilistic  one  and 
describe  the  cone  of  harmonic  functions  for  a  certain  diffusion.  We  adopt  the  notation 
from  [3]  as  is  customary  in  Markov  process  theory.  So,  the  state  space  is  still  locally 
compact  with  countable  basis,  but  now  it  is  denoted  by  E.  In  addition  we  assume 
that  E  is  connected.  Let  S  be  the  Borel  cr-algebra  on  E. 

We  assume  that  X  —  (£l,J-,J-t,Xt,6t,Px)  is  a  transient  diffusion  on  the  state 
space  (E,£).  By  this  we  mean  that  A"  is  a  transient  Hunt  process  with  continuous 
paths.  Let  T4  and  r^  denote  the  hitting  time  of  a  set  A,  and  the  exit  time  from  A, 
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respectively,  i.e., 

TA  =  inf{*  >  0  :  Xt  G  A}    and   ta  =  inf{*  >  0  :  Xt  <£  A}. 
Transience  means  that  for  every  compact  subset  K  of  E,  and  for  every  x, 

lim  PX{TK  o  6t)  =  0. 


t— >oo 


Let  (Pt)  denote  the  transition  semigroup  of  X  and  let  U  be  the  potential  oper- 
ator defined  by 

Uf{x)  =  /     Ptf(x)  dt 
Jo 

for  /  nonnegative  Borel  measurable  function.  Following  [8]  and  [21]  we  assume  that 
there  is  a  reference  measure  £  on  (E,  £)  such  that  the  potential  kernel  U(x,dy)  has 
the  density  u(x,  y)  with  respect  to  £.  That  is 

U{x,dy)  =  u{x,y)£(dy). 

It  is  assumed  that  u(x,  y)  satisfies  the  following  two  properties: 

(i)  (x,y)  t-+  u_1(x,j/)  is  finite  and  continuous,  and 
(ii)  u(x,  y)  =  oo  if  and  only  if  x  =  y. 

Note  that  (i)  implies  that  u(x,y)  >  0. 

A  nonnegative  measurable  function  /  is  excessive  if  (i)  Ptf  <  f  for  every  t  >  0, 
and  (ii)  limt_0  Ptf  =  /• 

For  a  Borel  set  A  and  a  measurable  function  /,  let 

PAf(x)  =  E*V(Xta);Ta  <  oo] 

be  the  hitting  operator  of  A.  A  nonnegative  measurable  function  h  is  called  harmonic 
(for  X)  in  an  open  subset  D  of  E,  if 

h(x)  =  PKch(x)  (3.18) 
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for  every  compact  subset  K  of  D.  If  D  =  E,  then  we  say  that  h  is  harmonic  in  E. 

It  was  proved  in  [8,  p. 127,  Thm.6],  that  every  excessive  function  s  which  is  not 
identically  infinite  has  the  Riesz  decomposition 

s  =  h+      u(x,y)fi(dy)  (3.19) 

Je 

where  h  is  excessive  and  harmonic,  and  ft  is  a  Radon  measure  on  (E,  £).  Furthermore, 
s  is  harmonic  in  D  if  and  only  if  fi(D)  =  0  (cf.  [21,  p.610,  Thm.2]. 

The  result  which  makes  harmonic  function  tractable  in  view  of  Chapter  2  is  the 
following  one  [21,  p. 610,  Thm.3]. 

THEOREM  3.1  Let  h  be  an  excessive  function,  h  ^  oo.  If  h  is  harmonic  in  an  open 
set  D,  then  h  is  continuous  in  D. 

We  outline  the  proof  following  [21]. 

PROOF:  We  assume  that  D  is  relatively  compact  open.  The  function  s  =  P^h  is 
excessive.  Since  5  =  h  in  D,  the  continuity  of  paths  of  X  implies  that  s  is  harmonic 
in  D.  Let  s  =  Ufi:  then  (i  is  carried  by  D  and  since  s  is  harmonic  in  D,  fi(D)  =  0. 
Further, 

s(x)  =  PDh(x)=   I    u(x,y)n(dy).  (3.20) 

JdD 

Continuity  assumption  for  u  yields  that  the  right-hand  side  is  continuous.  Hence,  h 
is  continuous  in  D.O 

With  the  same  assumptions  as  in  [21],  it  is  proved  in  [22,  Cor.l],  that  every 
nonnegative  harmonic  function  for  X  is  automatically  excessive.  So,  the  theorem 
above  says  that  all  nonnegative  harmonic  functions  are  continuous.  Let  7{  denote 
the  family  of  all  nonnegative  functions  which  are  harmonic  for  X  on  the  whole  space 
E.   Then  H  is  a  convex  cone  of  nonnegative  continuous  functions.    Every  function 
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in  Ji  is  strictly  positive  (unless  identically  zero).  Indeed,  let  h  €  W  and  h(xQ)  =  0 
for  some  Xo  E  E.  Let  D  be  a  relatively  compact  open  subset  of  E  containing  x0. 
If  s  =  Poh,  then  s  —  fdD  w(x,  y)n(dy)  for  a  Radon  measure  fi  concentrated  on  dD. 
Further,  s  =  h  in  D,  hence  0  =  s(x0)  =  JdDu(x0,y)fi(dy).  But  u(x0,y)  >  0  implies 
that  the  measure  \i  is  zero.  Therefore  5  =  0  and,  in  particular,  h  —  0  in  D.  Since  E 
can  be  covered  by  a  sequence  of  relative  compact  open  subsets  of  E,  it  follows  that 
h  -  0  in  E. 

Following  the  pattern  from  sections  3.1  and  3.2,  we  prove  that  %  is  closed  in 
C(E). 

Proposition  3.1   The  cone  H  is  closed  in  C(E). 

PROOF:  Let  {hn}  be  a  sequence  of  functions  in  Ti  which  converges  to  a  function  h 
uniformly  on  compacts.  Let  K  be  a  compact  subset  of  E.  If  x  £  K  is  regular  for 
Kc,  i.e.,  if  Px(TKc  =  0)  =  1,  then  XTh.c  =  x  Px&.e.,  so  PK'h(x)  =  h(x).  Note  that 
since  Kc  is  open,  every  point  of  Kc  is  regular  for  Kc.  Therefore  irregular  points  for 
Kc  are  only  in  K. 

Let  x  be  irregular  for  Kc,  i.e.,  Tk*  >  0  Pxa,.e.  Then  x  G  K  and  by  continuity 
of  X  we  have  XjKC  6  OK,  Pxa..e.  For  e  >  0,  we  can  find  n0  large  enough  such  that 
\hn(y)  —  h(y)\  <  e  for  all  y  €  K  and  all  n  >  n0.  Further, 

\PKchn(x)  -  PKch(x)\  <  Ex[\hn(XTKC)  -  K{XTkc)\-TKc  <  oo]  <  ePx[TKc  <  oo]  <  e, 

for  n  >  no-  Therefore, 

\h(x)  -  PKch(x)\  <  \h(x)  -  hn(x)\  +  \hn(x)  -  PKchn(x)\  +  \PKchn(x)  -  PK<h(x)\  <  2e, 

since  the  middle  term  is  zero  by  assumption.  This  proves  that  h  is  harmonica 
Now  we  fix  a  point  x0  in  E  and  consider  the  family 

Hxo  =  {k:h€  H,  h(x0)  =  1}.  (3.21) 


42 


Then  T-Cxo  is  obviously  a  basis  for  the  cone  7i.  Our  next  goal  is  to  show  that  it 
is  compact  in  C(E).  This  follows  from  Theorems  4  and  5  in  [21].  Theorem  4  shows 
that  Harnack  inequality  holds,  while  Theorem  5  gives  the  local  equicontinuity  of  7iXo. 
We  take  another  approach  and  prove  directly  that  7iX0  is  compact  without  invoking 
Harnack  inequality. 

Proposition  3.2  HXo  is  compact  in  C(E). 

PROOF:  Let  {hn}  be  a  sequence  in  ?^Xo.  We  want  to  show  that  there  is  a  subsequence 
{hnk}  and  h  in  7iXo  such  that  hnk  —*  h'm  C(E). 

Let  K  be  a  compact  subset  of  E  and,  without  loss  of  generality,  we  assume  that 
xo  6  K.  Let  D  be  a  relatively  compact  open  set  in  E  containing  K.  For  each  n  G  N 
we  define  vn  =  PDhn.  Then  by  (3.20) 


vn(x)  =  I     u(x,y)fin(dy) 

JdD 


where  fin  is  a  positive  measure  on  dD.  Further,  Pohn  —  hn  in  D. 

Since  y  i— >  u(x0,  y)  is  continuous  and  positive  on  dD,  it  attains  its  minimum, 
say  m  >  0.  Then 

1  =  hn(x0)  =  PDhn(xo)  =  /     u(x0,y)fin(dy)  >   /     m/jJdy). 

JdD  JdD 

Hence,  fJ.n(dy)  <  1/m  for  every  n,  and  thus  {fin}  is  a  bounded  sequence  of  measures. 
Let  {fink }  be  its  weakly  convergent  subsequence  and  fi  the  limit  measure,  which  is 
also  concentrated  on  dD. 

Let  v  =  Ufi  =  fdDu(-,y)n(dy).  Then  v  is  excessive  in  E  and  harmonic  in  D 
(since  fi(D)  -  0).  For  each  fixed  x  G  D,  the  function  y  »-*•  u(x,y)  is  continuous  on 
dD  and  therefore 

I     u(x,y)fink(dy)  ^   I    u(x,y)[i(dy), 

•JdD  JdD 
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i.e.,  hnk(x)  — »  v(x),  for  every  x  €  D. 

Now  we  show  that  this  convergence  is  uniform  on  K.  The  function  (x,  y)  \-* 
u(x,  y)  is  uniformly  continuous  on  K  x  dD.  Hence,  for  e  >  0,  there  is  8  >  0  such  that 
\u(x,  y)  —  u(x',  y')\  <  e  whenever  d(x,  x')  +  d(y,  y')  <  8  ( d  is  a  metric  compatible  with 
the  topology  on  E).  By  compactness  of  K,  there  are  finitely  many  balls  B(xi,8), 
i  =  1, . . . ,  /,  covering  K.  Let  k0  E  N  such  that  for  k  >  A;0,  \hnk  (x,)  —  u(xt)|  <  e  for 
i  =  1, . . . ,/.  For  a  given  x  £  if  we  choose  the  ball  B(xi,8)  containing  x.  Then,  for 
k  >  fco, 

<  \Unnk(x)  -  Unnk{xi)\  +  \Unnk(xi)  -  Un(xi)\  +  \Ufi(xi)  -  Un(x)\ 
<         \u{x,y)-u(xi,y)\nnk(dy)  +   \hnk(Xi)-v(xi)\   +    f    \u(xi,y)  -  u(x,y)\fi(dy) 

JdD  JdD 

<  efink(dD)  +  e  +  tfi{dD)  <  e(2/m  +  1). 

Therefore,  the  convergence  is  uniform  on  K. 

So  far  it  is  proved  that  for  every  sequence  {hn}  in  Wxo,  for  every  compact  set  K 
in  E  and  every  relatively  compact  open  set  D  containing  K,  there  exist  a  subsequence 
{hnk}  and  an  excessive  function  v,  harmonic  in  D,  such  that  hnk  -*  v  uniformly  on 
K. 

Let  {Kn}  be  an  increasing  sequence  of  compact  sets  covering  E  such  that  Kn  C 
IntA'n+1  and  x0  G  K1.  Let  {hn}  be  a  sequence  in  Hxo.  For  K\  there  is  a  subsequence 
{hnl}  of  {hn}  and  an  excessive  function  Vi  harmonic  in  Int/^2,  such  that  /inl  -^  ua 
uniformly  on  Kx.  We  proceed  inductively:  For  jeN  there  is  a  subsequence  {hnJ+1} 
°f  {hn,j}  and  an  excessive  function  vj+1,  harmonic  in  Int/^+2  such  that  hnj+1  -^  vj+1 
uniformly  on  Kj+\.  Obviously,  vj+i\Kj  =  Vj. 

Let  {hn<n}  be  the  diagonal  subsequence  and  define  the  function  v  on  E  by 
v\Kj  =  Vj.  Then  hn,n  ->  v  uniformly  on  Kj  for  every  j  e  N.  This  implies  that  v  is 


4-1 


continuous  and  that  hn<n  -*  v  in  C(E).  Since  every  function  hn,n  is  in  W^,  the  proof 
of  Proposition  3.1  implies  that  v  G  "H.  Trivially,  v(x0)  =  1.  □ 

Together,  the  last  two  propositions  imply  that  7{  is  a  closed  convex  cone  with 
compact  basis  'Hxo.  Functions  in  H  are  strictly  positive  (unless  zero).  By  Corollary 
2.1,  we  see  that  functions  in  7{  satisfy  Harnack  inequality.  The  family  7ilni  of  all 
possible  infima  of  functions  in  H  consists  of  strictly  positive  functions.  Every  function 
in  7iloS  is  excessive.  Indeed,  let  v  =  inf  a  hQ  be  a  function  in  Hinf  and  let  D  be  a 
relatively  compact  open  subset  of  E.  Then,  for  every  a  and  every  x  6  E, 

Ex[v(XTD)\  <  Ex[hQ(XTD)}  =  ha(x) 

which  implies 

E*[v(XTD)]<v(x). 

Since  v  is  continuous,  Theorem  5.11(11)  from  [3]  implies  that  v  is  excessive. 

By  the  Riesz  representation  theorem,  (3.19),  each  v  €  7iin£  can  be  written  as 
v  =  h  +  Ufi,  where  h  is  harmonic  and  (i  is  a  Radon  measure.  Obviously,  v  —  h€  "Hinf , 
and  hence  Ufi  €  7iinS. 

For  the  next  proposition  we  need  to  assume  that  the  state  space  E  satisfies  a 
certain  topological  property.  We  note  that,  for  example,  star-shaped  domains  in  Rn 
have  this  property. 

PROPOSITION  3.3  Let  E  satisfy  the  following  condition:  For  every  compact  subset 
K  of  E  there  exists  a  relatively  compact  open  set  D  containing  K  such  that  E\D  is 
connected.  Let  v  =  Ufi,  v  ^  0,  be  a  potential  from  Hinf .  Then  fi  cannot  have  compact 
support. 

PROOF:  Assume  that  Supp^  is  compact.  Let  D  be  a  relatively  compact  open  set 
containing  Supp/i.  Then  /j,(E\D)  -  0  ,  so  v  is  harmonic  in  E\D.  Let  x0  be  any  point 
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from  E\D.  By  Lemma  2.7,  there  is  h  £  W  such  that  h  >  v  in  E  and  /j(x0)  =  v(x0). 
Then  the  function  h  —  v  is  nonnegative  and  harmonic  in  E\D  and  (/i  —  u)(:ro)  =  0. 
Since  E  \  D  is  connected,  the  same  argument  as  before  shows  that  h  —  v  must  be 
identically  zero  on  E  \  D,  i.e.,  h  =  v  on  E  \  D.  By  continuity,  h  =  v  on  3D.  Then 
we  have,  for  any  x  in  D, 

h(x)  =  E*[h{XTD)\  =  E*[v{XTD)\  <  v(x) 

since  XTD  £  3D.  Therefore,  h  <  v  in  D.  But  h  >  v  everywhere  so,  h  —  v  in  E.  This 
is  a  contradiction,  because  a  potential  cannot  be  harmonic.  □ 

If  the  assumption  on  E  is  not  satisfied,  then  it  can  happen  that  fi  has  compact 
support.  To  see  this,  let  E  be  the  unit  disc  in  R2  punctured  at  the  origin:  E  =  {x  € 
R2  :  0  <  |x|  <  1}.  Let  u  =  (—  log  \x\)  A  1.  Since  —  log  \x\  is  harmonic  and  positive  in 
E,  u  G  7Yinf(£).  But  the  Riesz  measure  of  u  sits  on  the  circle  centered  at  the  origin 
with  radius  1/e,  so  it  thus  has  a  compact  support. 


CHAPTER  4 
INFIMA  OF  HARMONIC  FUNCTIONS 


4.1     Function  k  for  the  Unit  Ball 

In  this  chapter  we  consider  in  more  detail  the  family  of  all  infima  of  nonnegative 
harmonic  functions  in  some  domain  D  in  R".  Here  harmonic  is  understood  in  the 
classical  sense:  h  :  D  — >  R  is  harmonic  if  Ah  —  0  where  A  denotes  the  Laplacian. 
Let  Ti  denote  the  family  of  all  nonnegative  harmonic  functions  in  D,  and  W  is  as 
usual.  Then  the  discussion  in  3.1  shows  that  H.  is  a  closed  convex  cone  with  compact 
basis,  and  when  D  is  a  bounded  Lipschitz  domain  there  is  a  kernel  function  satisfying 
(Hi)-(H3)  and  so  the  results  from  Chapter  2  are  applicable. 

In  order  to  compute  something,  we  restrict  ourselves  to  a  very  special  domain 
—  the  unit  ball  in  Rn.  By  using  the  explicit  form  of  the  Poisson  kernel  and  the 
Mobius  transform,  it  is  possible  to  get  some  explicit  formulae. 

Let  D  be  the  unit  ball  in  Rn.  Recall  the  definition  of  functions  sa,  a  G  D: 

sa(x)  =  k(x,a)  =  inf{^(x)  :  h  E  H,  h(a)  =  1}. 

First  we  compute  s0  where  0  denotes  the  origin. 

Every  positive  harmonic  function  in  D  satisfies  the  classical  Harnack  inequality 
(e.g.  [14,  p.31]) 


(T^(O)  <  Mx)  <  (rfjj*(0).  (4.D 


If  h(0)  =  1,  then  the  left  inequality  reads 

1-lxl2 


(l  +  l*l)n 
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<  h(x). 
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Therefore, 

l-H2 


<  s0(x).  (4.2) 


(1  +  |*|)- 

Let  P(x,z),  x  G  D,  z  £  dD,  denote  the  Poisson  kernel  for  the  unit  ball: 

1  -  Ixl2 
P(s,*)=,       '      .  (4.3) 

For  iGD,  let  z  and  —  z  denote  points  from  dD  which  lie  on  the  diameter  through  x 
with  —z  closer  to  x.  Then  \z  —  x\  =  1  +  |x|,  and  so 

P(x,z)  = ^—. 

(i  +  kl)n 

Together  with  (4.2)  this  yields 

In  order  to  get  5a  for  arbitrary  a  £  Z),  we  need  some  results  on  the  Mobius 
transforms.  These  can  be  found  in  [1]  and  [17]. 

A  similarity  in  Rn  is  the  mapping  x  (-►  Mx  +  6,  where  b  6  Rn  and  M  is  a 
conformal  matrix,  i.e.,  M  =  A/f  with  A  >  0  and  K  €  O(n).  Here  0(n)  denotes  the 
orthogonal  group.  It  is  a  trivial  fact  that  all  similarities  preserve  harmonicity. 

Reflection  through  the  unit  sphere  is  the  map  J  :  x  y-*  x*  defined  by  x*  = 
x/\x\2.  Following  [1],  let  Q(x)  denote  the  matrix  with  entries  Q(x)ij  =  XiXj/\x\2. 
The  derivative  of  J  is 

A*)  =  pjiC  -  2<?(*))  (4.5) 

(see  [1,  p.l8,(14)]),  where  /  denotes  the  identity  matrix.  Furthermore,  (I -2Q)2  =  I, 
i.e.,  I-2Qe  0(n). 

The  full  Mobius  group  M(Rn)  is  the  group  generated  by  all  similarities  together 
with  J.    For  any  7  €  M(Rn),  we  denote  by  |7'(rc)|  the  positive  number  such  that 
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l'(x)/\Y(x)\  €  0(n).  It  follows  from  (4.5)  that  \J'{x)\  =  l/\x\2.  If  7  =  XK  +  b  is  a 
similarity,  then  |7'(x)|  =  A. 

Let  h  be  a  positive  harmonic  function.  The  Kelvin  transformation  of  h  is  the 
function 

*(*)  =  ppi*  (pji)  • 

Kelvin  transformation  preserves  positivity  and  harmonicity  (e.g.,  [23,  4.25]).    Note 
that  h  can  be  written  as 

h  =  IJ'ix^hiJx). 

Hence,  if  7  is  either  a  similarity  or  the  reflection  «/,  then  the  map 

T^-.h^h,        h  =  |7'(x)|IT2/i(7x)  (4.6) 

preserves  positivity  and  harmonicity.     By  the  chain  rule,  this  holds  for  any  7  6 
M(Rn).  Hence  the  following  lemma  from  [17]  holds: 

LEMMA  4.1  For  any  7  6  M(Rn),  the  map  T7  defined  by  (4.6)  preserves  positivity 
and  harmonicity. 

We  shall  need  a  Mobius  transformation  which  maps  the  unit  ball  onto  itself. 
For  a  =  0,  let  70  be  the  identity,  and  for  a  /  0, 

~  (x)  -  (l-|«|2)(x-a)-|x-«|2a 

where  ax  =  a1x1  +  a2x2-] \-anxn.  Then  70  G  M(Rn),  7a  maps  D  onto  D,  ja(a)  =  0 

and  7_0  07,,  =  id^  (for  proofs  see  [1]).  Furthermore, 

,      .   .,  \x  —  a\ 

Mx*  =  dcrir  (4-8) 


and 


1   //   \i  1  —  lal2 
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(cf.  [1,  pp.26,  27,  formulae  (30)  and  (32)]). 

Now  we  define,  for  every  a  €  D,  the  map  Ta  :  h  \— *  h,  by 


la -a*"1"2 


M*)  = ±^M7«(*)),    xeD  (4.10) 


j^r^pi^.(*)).    *ei> 


for  a  ^  0,  and  To  =  id.  By  using  (4.9),  this  can  be  rewritten  as 

By  Lemma  4.1,  we  see  that  Ta  preserves  harmonicity  and  positivity.  Moreover,  if 
h(0)  =  1,  then  h(a)  =  1.  If  for  b  G  D,  Hb  =  {h  :  h  €  H,  h(b)  =  1},  then  the 
discussion  above  shows  that  Ta  :  Tio  —>  W0. 

Now  we  show  that  Ta  is  onto.  Let  u  €  "Ha  and  define 

&w=M^pu(7-'(x))-  (4J1) 

Then  /&  is  positive  and  harmonic  in  Z),  and  ft(0)  =  u(7_a(0))  =  u(a)  =  1.  Further, 

Tah{x)  =  jfr^*(7.(*)) 

|a-a1"-2Kr2 
=  |x-a|n-2|7a(x)  +  «r-2U(7-a(7a(:C))) 

n-2 

|        u(x)  =  u(:r). 


\J:r  -  a||7a(x)  +  a*|/ 

The  last  equation  follows  from  an  easy  computation  showing  that  the  factor  before 
u{x)  is  identically  1. 

Now  we  are  ready  to  compute  sa. 

sa(x)  =  mi{u(x)  :ueHa}  =  mf{(Tah)(x)  :  h  €  H0} 

{\a  -  a*\n~2  ~\         \n  -  n*\n~2 
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Now  we  use  (4.4)  and  (4.8)  and  get 

|a_«T-2     l_|7a(l)p 


*a(x) 


\x  _  a*|»-2  (i  +  |7a(x)|)« 


|a_a»|n-2   1-  |a|^t|2  i_i2xn-2  |a|2|ar  —  a*|2  —  |ar  —  a|2 

■  |*  -  a*|-»  (i  +  1_J^)"       l        '   N        (|a||*  -  a'|  +  |x  -  a|)»' 

Since  |a|2|x  -  a*|2  -  |x  -  a\2  —  (1  -  |a|2)(l  -  |x|2),  we  finally  obtain 

,  <v»     n     1.12^-2    (l-!«l2)(l-kl2)        „,n  ai9x 

S°(X)  =  (1  -  W  }       (|«||*-a-|  +  |*-a|)-'    ^°-  (°2) 

It  is  obvious  from  this  formula  that  lim  sJx)  =  0  for  every  z  6  5Z).  Therefore, 

s0  is  a  potential.  The  Riesz  measure  of  sa  is  the  measure  with  density  —  Asa  (with 

respect  to  the  Lebesgue  measure  on  D).  We  will  spare  the  reader  of  several  pages  of 

computation  and  simply  write  down  the  formulae: 

-As°{x)=\*t+wSo(x)'  x^°  (4'13) 

and,  for  a  ^  0, 

a     /   ^  n(n  -  1)(1  -  \a\2)n  ,   ,  ,  ,        , 

"  **W  =  \a\\X-a*\\z-a\{\a\\X-J\  +  \X-a\y  «"«■    *  *  *'  ^ 

According  to  Proposition  2.4,  for  any  x  6  D  there  exists  a  point  z  =  z(0,  x)  on 
the  boundary  dD  such  that  s0(x)  =  P(x,z).  It  is  obvious  from  the  way  how  s0  was 
computed  that  z  =  -x/|x|   for  x  ^  0,  i.e.,  z  is  the  point  on  the  diameter  through  x. 

If  we  are  given  z  6  3D,  then  the  set  of  points  x  in  D  for  which  s0(x)  =  P(x,  z) 
is  the  segment  [0,  -z)  =  {Xz  :  -1  <  A  <  0}.  Let  us  denote  this  segment  by  L0iZ. 
Thus,  L0,z  =  {x  €  D  :  s0(x)  =  P(x,  z)}. 

There  is  nothing  special  about  the  point  0,  so  we  fix  a  point  a  G  D,  a  ^  0. 
Again  by  Proposition  2.4,  for  any  x  €  D  there  is  a  point  £  =  ((a,x)  €  dD  such 
that  sa(x)  =  P(x,Q/P(a,().  To  find  such  f  we  go  backwards  .  We  fix  (  €  dD  and 
determine  the  set  La<<  =  {x  e  D  :  sa(x)  =  P(x,C)/P(a,Q}. 
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Let  7o  be  as  in  (4.7).  We  obtained  that 

la-«T-2     ,     (   vx 
5a(x)  =  \x  -  a*\"-*S°Mx)) 

for  every  x  €  D.  Let  2  =  7o(C)  G  ^  and  choose  x  such  that  7a(x)  £  L0<z  =  (0,  — z). 
Then  s0(7a(z))  =  ^(7a(z),2),  and  so 

l«-«T-2pr     f   ,       ,,»       I«-«1-2       l-|7„(*)|a 
MXj-  |«-*|^™*''™»-  |x-a-M   l%(«)  -  7.(01"' 

By  formula  (15),  p.19,  in  [1],  |7.(x)  -  7«(C)I  =  Wa(x)\1/2\Ya(0\1/2\x  -  CI-  A  simple 
computation  using  the  above,  (4.8),  and  (4.9)  determines  sa(x)  —  P{x,  ()/P(a,  (). 

This  shows  that  La^  =  7~1(L0,2)  =  7-a(Lo,z),  where  z  =  7a(C)-  Conformal 
property  of  -fa  easily  implies  that  La^  is  an  arc  of  the  circle  through  a,  perpendicular 
to  the  sphere  dD.  Precisely,  if  £  is  given,  let  S  be  the  circle  through  a  perpendicular 
to  the  sphere  dD.  Then  S  D  dD  =  {£,  £1}.  The  arc  of  S  connecting  a  and  (1  is  La^. 

Therefore,  if  x  6  Z),  let  5  be  the  circle  through  a  and  x  perpendicular  to  dD. 
The  point  (  €  dZ)  such  that  sa(x)  =  P(x,()/P(a,Cl)  is  the  point  on  S  D  dD  such 
that  £  does  not  belong  to  the  arc  connecting  a  and  (.  Note  that  £  is  unique. 

Let  us  denote  sa(x)  by  &(x,  a)  as  in  Chapter  2.  For  n  =  2, 

*l:r'Gj-(|a||z-a*|  +  |z-<z|)2 

which  is  symmetric  in  x  and  a  due  to  the  fact  that  |a||x  —  a*\  =  \x\\a  —  x*\.  When 
n  >  3,  k  is  noi  symmetric. 

One  can  consider  quotients  of  the  form  sa(x)/sQ(x).  Then  straightforward  com- 
putation shows  that,  for  every  z  €  dD, 

x~*z  s0(x) 
which  was  the  motivation  for  Section  2.2. 
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It  is  tempting  to  compare  the  function  k(x,a)  with  the  familiar  Green  function 
G(x,a)  for  the  unit  ball  D.  Recall  the  formulae  for  the  Green  function: 

1  —  xa 


G(x.a)  =  log 

x  —  a 

1 


n  =  2  (4.15) 


G(x,a)  =  -l r-T-    r-r— r: 1— r,    n>3  (4.16) 

v   '    '       |x-a|n"2        \x\n~2  \x  -  a*\n~2  ' 

where  a  is  the  conjugate  of  a,  and  a*  =  Ja  (e.g.  [23,  6.37]). 

Green  function  is  infinite  on  the  diagonal  while  k  is  1  on  the  diagonal.  Therefore, 
it  doesn't  seem  that  G  and  k  can  be  compared  away  from  the  boundary  dD.  So,  let 
us  consider  both  function  when  one  of  the  variables  is  near  the  boundary. 

First  we  let  n  =  2.  Then 

k(0,a)       l-|q|        1  _L bgjg| 

G(0,a)  ~  l  +  |a|log|l/a|  "  1  +  \a\  \a\  -  l' 

Since  (\a\  —  l)/log  \a\  — *  1  as  a  — >  2  €  <9Z),  it  follows 

For  igj),  both  G(x,a)/G(0,a)  and  k(x,a)/k(0,a)  converge  to  P(x,z)  as  a  — ►  2. 

Hence,  from  (4.16)  it  follows 

Km  ^l      ^^ 
a->*Cj(x,a)       2 

For  n  >  3  the  situation  is  different.  First  we  consider  k(x,a)/G(x,a)  as  x  — ► 
0  €  d.D.  Straightforward  computation  using  (4.12)  and  (4.16)  gives 

hMJ'-H'r;,  2£M3.  (4.19) 

*~»G(x,a)       2"-1(n-2)'  v        ; 

This  shows  that  when  x  is  close  to  the  boundary,  the  infimum  sa(x)  is  comparable 

with  the  Green  function  G(x,a).  On  the  other  hand,  since 

,.      k(x,a)       ,.      (l-|a|2)"-2       n 
hm  tJ-1-(  =  hm  )- '    n  .  =  0, 

'•-'^  fc(a,  x)       a~>2:  (1  -  |x|2)"-2 
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it  follows  that 

k(x, a)       ,.      k(x,a)  k(a,x)       „  ,._  ,.  „„» 

lim-^-+  =  lim  V       ;=0,    zedD.  4.20 

a—zG(x,  a)       °--J:  k(a,x)  G(a,x) 

Therefore,  for  a  fixed  x  G  Z),  when  a  approaches  the  boundary,  s„(x)  goes  to  zero 
faster  then  G(x,a).  This  suggests  that  the  fact  that  k(x,a)/k(x0,a)  tends  to  the 
kernel  function  K(x,  z)  as  a  — ►  z  G  dD  (here  Z)  is  Lipschitz,  and  xQ  G  Z)),  is  not 
obtainable  from  the  corresponding  result  for  the  Green  function. 

We  also  note  that  k(x,a)  ~  6(a)n~1  when  a  is  near  the  boundary  dD,  and 
k(x,a)  ~  8(x)  when  x  is  near  5.D,  where  <$(•)  denotes  the  distance  to  the  boundary. 

4.2     Function  fc  for  Some  Other  Domains 

In  this  section  we  use  transformations  that  preserve  harmonicity  to  obtain  ex- 
plicit formulae  and  some  properties  of  sa  for  certain  domains. 

Let  B(x0,  r)  be  the  ball  in  Rn  centered  at  x0  with  radius  r.  The  reflection 
through  the  sphere  S(x0,r)  =  {x  G  Rn  :  \x  —  xQ\  =  r}  is  the  mapping  J  defined  by 

r2 

J(x)  =  x*  =  x0  +  -. ^(x-x0).  (4.21) 

|x-x0|2 

Obviously  J2  —  idj^". 

Let  D  be  an  open  set  in  Rn  \  {x0}  and  let  D*  =  J(D).  The  Kelvin  transform 

of  a  function  h  on  D  is  a  function  h*  on  D*  defined  by 

ft"(*-)  =  ]~^  *(*)■  (4.22) 

Kelvin  transformation  preserves  positivity  and  harmonicity  (e.g.,  [23,  4.25]).  Let 
H(D)  and  H{D")  denote  families  of  positive  harmonic  functions  on  D  and  D*  re- 
spectively. For  a  G  D  we  define  the  mapping  Sa  :  H(D)  -»  H(Dm)  by 


{s"k)ix') = S-ZF k(x)-  <4-23) 
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Then  Sah  6  7i(D*)  and  (Sah)(a*)  =  1.   It  is  easy  to  see  that  5„  is  1  —  1  and  onto. 
Therefore,  as  in  the  previous  section,  we  get  the  formula 

Further,  in  the  same  way  as  in  Section  4.1,  one  can  compute  the  function  sa  for 
the  ball  B(a,p).  It  follows 

*(•>  =  '"-2fe^p  <4-25> 

Now  we  compute  the  function  sa  for  the  upper  halfspace.  We  work  in  the 
(n+  l)-dimensional  Euclidean  space  Rn+1.  Points  in  Rn+1  are  denoted  by  P,  Q,R,.. ., 
P  =  (x,  y),  x  6  Rn,  y  6  R.  Let  #+  =  {(x,  t/)  6  Rn+1  :  y  >  0}  be  the  upper  halfspace. 
We  fix  the  point  Q*  G  H+,  Q*  =  (0, v),  v  >  0.  The  goal  is  to  compute  sq>. 

Let  «7  denote  the  reflection  through  the  sphere  S(—Q*,v),  and  let  D  be  the  ball 
B(—Q*/2,v/2).  Then  D*  =  H+  and  an  easy  computation  gives  J(—Q*/2)  =  Q*. 
Let  P*  g  £f+;  then,  by  (4.24)  and  (4.25), 


2|Q*|        /^V"1     I  —  Q*/2|2  —  |,P  H-  Q*/2| 


2 


\P  +  Q*!""1  V27       (|  -  Q-/2|  +  |P  +  Q*/2|)n+1 ' 
After  an  easy  computation  one  gets 

,     fm_         ^"~2  ^/4-|P  +  Q72|2 

on  ;  -  |P.  +  q.|«-i  (r/2  +  |p  +  g*/2|)n+1 '  l     J 

Let  P  =  (x,y),   x  6  Rn,  y  >  0.    Straightforward  computations  yield  the  following 
formulae: 

|P*  +  Q*|  =  [H2  +  (2,  +  r,)2]1/2  (4.27) 
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\p+<rm'=T^Trr^  <4-29> 

4  \x\2  +  (y  +  v)2 


l^  +  ^/2|2=,    ,"",.,„.  (4-30) 


By  putting  together  (4.27)-(4.30)  -  and  after  some  simplification  -  we  finally  obtain 


sQ.(P*)  =      ,  = ,  y (4.31) 

(J\x\2  +  (y  +  v)2  +  y/\x\*  +  (y  -  v)*)»» 

where  (P*)  =  (x,y)  and  Q*  =  (0,v). 

If  Q*  =  (u,v),  u  £  Rn,  then  one  should  replace  x  by  x  —  u  m  (4.31)  to  obtain 
the  formula  for  Sq*.  When  we  are  in  the  two-dimensional  space,  that  formula  gives 
the  symmetry  in  P*  and  Q*. 

One  can  again  consider  the  limit  of  the  quotient  sq*(P*)/sq.(P0*)  as  Q*  tends 
to  some  point  (2,0)  on  the  boundary  of  H+.  Here  P0*  =  (x0,y0)  is  some  fixed  point 
in  H+.  The  limit  is  equal  to 

y(\x0-t\2  +  y02Yn+^2 
2/o  (Ix-^  +  y2)^1)/2  ' 

which  is  a  constant  multiple  of  the  kernel  function 

for  the  upper  halfspace  H+ .  Finally  we  note  that  as  Q*  goes  to  infinity  (i.e.,  v  — >  oo), 
the  quotient  sq*(P*)/sq*(P0*)  approaches  y/yo. 

Now  we  consider  domains  in  R2.  We  identify  R2  with  the  complex  plane  C. 

Let  D  =  {z  G  C  :  \z\  <  1}  be  the  unit  disc  in  C  and  let  Di  be  a  simply 
connected  region  in  C  (with  at  least  two  boundary  points).  Then  there  is  a  conformal 
mapping  w  :  D\  — *  D.  Let  w  =  u  +  iv.  For  /  :  D  — *  R  of  class  C2,  we  define 
flf :  2>i  -♦  R  by 

g(x)  =  f(w(x)),     x  G  Dx.  (4.32) 
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Then  g  G  C2(DX)  and 

Ag(x)  =  \Vu(x)\2(Af)(w(x)).  (4.33) 

(e.g.  [23,  6.19]). 

If  /  is  positive  and  harmonic  in  D,  i.e.,  /  G  T-C(D),  then  g  is  positive  and 
harmonic  in  D\,  g  G  H(Di).  This  follows  from  definition  and  (4.33).  Moreover, 
g(x)  =  1  if  and  only  if  f(w(x))  =  1. 

For  a  G  Du  let  5.(y)  =  inf{flf(y)  :  y  G  W«(Z>i)}.  Then 

«.(*)  =  M{f(w(x))  :  /  G  Hw(a)(D)}  =  «^.)(w(«)).  (4.34) 

Therefore,  once  w  is  known,  one  can  obtain  the  formula  for  sa. 

Let  k(x,a)  =  sa(x).  Then,  from  (4.34),  we  see  that  k  is  symmetric  in  x  and  a. 

4.3     Examples 

It  is  not  easy  to  find  nontrivial  examples  of  functions  in  ?i,nf  on  a  given  domain 
D.  By  trivial  we  mean  functions  that  are  explicitly  written  as  an  infimum  of  positive 
harmonic  functions.  The  other  class  of  functions  in  Hini  which  may  be  regarded  as 
trivial  are  positive  concave  functions.  Yet,  they  are  useful,  because  we  can  use  them 
to  get  some  nontrivial  examples. 

Note,  on  the  other  hand,  that  we  do  have  a  necessary  condition  for  a  potential 

to  be  in  Ti"*1.  Let  Gd  denote  the  Green  function  on  D,  and,  for  a  positive  measure 

\i  on  D,  let 

GDn(x)=  /  GD(x,y)n(dy). 
Jd 

If  Gdh  G  H™*,  then  \i  cannot  have  compact  support  (Proposition  3.3).  Of  course, 
this  condition  is  far  from  being  sufficient.  There  are  measures  without  compact 
support  whose  potentials  are  not  even  continuous.  It  seems  that  the  characterization 
of  measures  which  give  potentials  in  H"**  is  difficult,  if  not  impossible. 
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In  this  section  we  give  some  examples  of  functions  in  riiat.  We  start  with  an 
auxiliary  result. 

Let  Dx  C  Rni,  D2  C  Rn2,...,  Dk  C  Rnfc  be  bounded  domains  and  let 
D  =  Dx  x  D2  x  . . .  x  Dk  be  a  bounded  domain  in  R"i+"2+-n*.  Assume  that 
U{  G  Winf(Z),),  i  =  1,2, . . . ,  k.  Then  the  function  u  defined  by 

u(x)  =  u(xi,  x2,  ...,«*)  =  ui(a;i)tt2(a;2)  •  •  •  ujt(zfc) 

is  in  riint(D).  To  see  this,  we  fix  x  in  ZJ.  For  each  i,  there  is  a  harmonic  function  h{ 
on  Di  such  that  /i,  >  u,  and  /ii(xj)  =  u,(a:,)  (see  Lemma  2.7).  We  define  h  on  Z)  by 

MS/)  =  h(yuV2,---,yk)  =  h1(y1)h2(y2)...hk(yk). 

Then  /&  is  harmonic  on  D,  h  >  u,  and  /i(a;)  =  u(a;).  Since  such  an  h  can  be  found  for 
any  x  in  D,  it  follows  that  u  6  HinI(D). 

The  above  can  be  used  to  get  some  examples  on  product  domains  in  Rn.  We 
assume  that  D  =  («j,  5i)  X  •  •  •  X  (an,  6n).  Let 


&(j,-)  =  sin  hr   '        ' j  ,  t'  =  l,...,n. 


Then  fc  is  the  eigenfunction  for  d2/dxf  corresponding  to  the  eigenvalue  — 7r2/(6,— a,)2. 
Since  <fn  is  concave,  it  is  in  ?iinf(a,,  6,).  Let  <f>  be  defined  on  D  by 

<£(:r)  =  </>(xx,  ...,xn)-  <f>i(xi) . . .  <f>n(xn). 

Then  (j>  is  the  eigenfunction  of  A  corresponding  to  the  eigenvalue  Y^=\{~ k2 l(bi— a,)2). 
Since  <f>  is  positive,  it  is  the  first  eigenfunction  of  A  on  D.  By  the  result  above,  (f>  is 
in  7{,ni(D).  Note  that  <f>  is  not  a  concave  function. 

Let  (Xt,  Px)  be  the  Brownian  motion  on  Rn  and  let  to  denote  the  first  exit 
time  from  D, 

td  =  mt~{t  >0:XtgD}. 
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Let  tj)(x)  =  Px(td  >  t).  It  is  well-known  that  0  is  excessive.  We  show  that  tp  is  in 
Winf(Z)). 

Let  r,  =  mi{t  >  0  :  Xj  £  (a,,  6,)},  i  =  l,...,n,  where  Xt  denotes  the  i-th 
component  of  A".  Since  X^  are  independent  processes,  times  t,-  are  independent. 
Further,  to  =  T\  A  . . .  A  rn.  Hence, 

0(x)  =  P*(rD  >  t)  =  P*(r,  >  0 .  ..P*(rn  >  0  =  ^(^i  >  0  •  •  -^"K  >  *), 

where  x  =  (xi, . . .  ,  x„).  But  each  X{  i-»  Px'(ri  >  t)  is  excessive  for  the  1-dimensional 
Brownian  motion  X^,  and  is  therefore  an  element  of  H"1* (a,-,  6,-).  Since  tj;  is  a  product 
of  such  functions,  it  follows  that  xj>  is  in  W    (D). 

Let  Z)  be  any  bounded  domain  in  Rn  and  let  rD  be  the  exit  time  from  D  of  the 
n-dimensional  Brownian  motion.  We  show  that  the  function  <f>(x)  =  Ex(t£>),  i.e.,  the 
expected  exit  time  from  D  is  in  7iinf(Z)). 

First  we  establish  this  result  for  the  ball  B  =  B(Q, r).    Let  tb  =  inf{i  >  0  : 

Xt  <?  B}.  Then 

Ex(tb)  =  -(r2  -  |z|2),     xeB 
n 

(e.g.  [23,  4.6]).  The  function  x  i-»  r2  —  |x|2  is  concave  on  Z?,  and  therefore  in  Hini{B). 
Now  we  take  up  the  general  case.  Since  D  is  bounded,  there  is  a  ball  B  =  B(0,  r) 
such  that  D  C  B.  Let  <f>{x)  =  Ex(td),  x  G  D,  and  V(z)  =  Ex(tb),  x  G  £.  It  is 
known  that  A</>  =  —1  in  D,  and  At/>  =  —1  in  B.  Hence  A(t/>  —  <£)  =  0  in  D,  so  0  —  <j> 
is  harmonic  in  D.  Obviously,  cf>  <  i\>  in  D.  If  h  =  ip  —  (f>,  then  /i  is  strictly  positive 
and  harmonic.  Since  ip  is  in  Ti^lB),  certainly  tp  G  7i,nf(Z)).  Hence,  ^  =  inf Q  /ia,  ha 
harmonic  and  positive  in  D.  But  then  <f>  =  i\>  —  h  —  mfa(ha  —  h)  and  each  ha  —  h  is 
positive  and  harmonic  in  D. 
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Let  (Pt)  be  the  semigroup  of  the  Brownian  motion  in  Rn.  Then,  for  nonnegative 
measurable  /, 

Ptf{x)  =  JRnp(t,x,y)f(y)dy, 

where 

P(^'y)  =  (^exp(-lfir-)- 

It  is  easy  to  see  that  the  function  y  i— ►  exp(— \y\2/2t)  is  concave  in  the  ball  5(0,  y/i). 
Therefore,  y  t->  p(t,x,y)  is  concave  in  B(x,\/i)  and  x  h-»  p(t,x,y)  is  concave  in 
B(y,Vt). 

Let  us  fix  T  >  0  and  let  t  >  AT2.  U  y  €  5(0,  T),  then  |x  -  y\  <  2T  for  each 
x  G  5(0,  T).  Therefore  x  G  5(y,2T)  C  B(y,\/i).  This  proves  that  x  •->  p(t,x,y)  is 
concave  in  5(0,  T)  for  each  ?/  G  5(0,  T).  In  particular,  a:  i-»  p(t,x,y)  is  an  element 
of  7iinf  on  5(0,  T)  for  any  y  G  5(0,  T). 

Let  p.  be  a  positive  Radon  measure  on  5(0,2).  By  using  a  result  from  2.1  we 
see  that  the  function  x  i-»  fB(o,T)P(^ix^y)^(^y)  1S  m  ^"^  on  B(®iT).  In  particular, 
if  /  is  nonnegative  and  /  G  £}OC(5(0,T)),  then  x  •-►  fB(0<T)p(t,x,y)f(y)(dy)  G 
ftinf(5(0,T)). 

Let  D  be  a  bounded  domain  in  Rn  and  let  /  be  a  nonnegative  function  which 
is  identically  zero  outside  D  and  /  G  C\  (D).  Let  T  be  large  enough  so  that 
D  C  5(0,  T)  and  let  t  >  AT2.  Then,  for  x  G  5, 

ptf(x)=  LnP(i,x,y)f(y)dy=   I         p(t,x,y)f(y)dy, 
JK  Jb(o,t) 

and  from  above  it  follows  that  Ptf  G  HinS(D). 
Further,  t  i— >  Ptf(x)  is  continuous,  so 


./4T2 


f// 
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is  in  7i,nf(Z))  (provided  that  it  is  not  identically  +00).  Let 

Gf(x)  =  /     Ptf{x)  dt. 
Jo 

Then 

I"0  Ptf(x)  =  P47^Gf(x)  =  GP4I*f(x). 

J4T2 

Hence  GPAT*f  €  Ww(i)),  and  similarly,  for  each  t  >  4T2,  GPtf(x)  G  ftinf(£). 

Let  Gd  denote  the  Green  function  for  D.  If  u  is  nonnegative  function,  then  it  is 
well  known  that  Go  u  =  Gu  —  h,  where  h  is  a  positive  harmonic  function  on  D  (e.g. 
[23,  6.2]).  By  using  this  fact  and  the  above  result,  it  follows  that  for  nonnegative 
/  e  C\oc(D)  there  exists  T  >  0  such  that  for  every  t  >  AT2,  GD(Ptf)  G  Hild{D). 

4.4     Two  Counterexamples 

In  this  section  we  give  two  examples.  One  of  them  shows  that  7i,nf  is  not  a 
potential  cone,  and  the  other  that  Winf  is  not  a  function  cone.  In  both  cases  we  take 
for  our  domain  the  unit  disc  in  R2. 

To  say  that  "Hinf  is  not  a  potential  cone  is  rather  vague,  since  we  have  not 
defined  what  is  meant  by  this.  There  are  different  approaches  and  axiomatizations 
of  potential  theory.  We  shall  adopt  the  one  from  [2]  and  show  that  H,,nt  is  not  a 
balayage  space.  Following  [2,  pp.  56-57],  we  consider  a  domain  D  and  a  convex  cone 
W  of  positive  lower  semicontinuous  functions  on  D.  The  (W— )  fine  topology  on  D 
is  the  coarsest  topology  on  D  which  is  finer  than  the  initial  topology  and  for  which 
all  functions  in  W  are  continuous.  Then  (D,  W)  is  said  to  be  a  balayage  space  if  the 
following  axioms  are  satisfied: 

(Bi)  If  {un}  is  an  increasing  sequence  of  functions  in  W,  then  limnun  G  W. 

(B2)   inf  ua    G  W  for  any  subfamily  {uQ  :  a  e  A},  where  "*  denotes  the  fine  l.s.c. 

regularization. 


61 

(B3)  K  u,Vi,v2  £  VV  such  that  u  <  v\  +  v2  there  exist  Ui,u2  €  W  such  that 

u  =  Ui  -f  «25  «i  <  ui5  «2  <  v2. 

(B4)  There  exists  a  function  cone  V  C  C+(Z>)    such  that  every  function  u  G  VV  is  an 

increasing  limit  of  functions  in  V. 

For  u  in  VV  and  any  subset  A  of  D  let 

iEjJ  =  inf{u  e\V  :v>uon  A}. 

It  is  proved  in  [2,  p.243,  Prop. 1.1],  that  ii  A  C  D  and  u,v  G  VV,  then 

*Sn  =  #  +  #  •  (4-35) 

We  take  for  D  the  unit  disc  in  R2  and  consider  the  convex  cone  Tiini  on  D. 
Then  (D^H^)  satisfies  axioms  (B1),(B2)  and  (B4)  (for  (B4)  see  Lemma  2.9).  Since 
(B3)  is  difficult  to  disprove  directly,  we  show  that  (4.35)  does  not  hold. 

Let  a  =  (-1/2,0),  b  =  (1/2,0)  be  two  points  in  D  and  let  A  =  {a, b}.  Then 
#£  =  inf{u  G  Ww  :  u(a)  >  sa(a),u(b)  >  sa(b)}  =  sa,  and  similarly,  R*  =  sb. 
Hence,  Rfa  +  Rfb  =  sa  +  sb.  Further, 

Rfa+sb  =  inf{u  E  H  :  u(a)  >  sa(a)  +  sb(a),u(b)  >  sa(b)  +  sb(b)} 

=  inf{u  €  H  :  u(a)  >  7,u(&)  >  7}, 

where  7  =  sa(a)  +  s&(a)  =  sa(6)  +  ^^(t)  because  of  symmetry.   Now  we  use  formula 
(4.12)  to  compute  sa  and  sb.  It  follows  that  sa(b)  =  1/9,  so  7  =  10/9. 
Let  x0  =  (0, 1/2);  then,  by  (4.12)  and  by  symmetry, 

sa(xo)  =  sb(x0)  =    ^-^     -      w  0.1862436,  (4.36) 


and  so 


R?a(x0)  +  <(x0)  =  aa(*0)  +  sb(xQ)  =  ^  w  ( 
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Now  we  find  a  function  h  G  H  such  that  h(a)  >  7,  h(b)  >  7,  but  h(x0)  < 
sa(xo)  +  sb(x0).  This  will  show  that  Rfa+,b(xo)  <  Rt(xo)  +  Rfb(x°)^  which  contradicts 
(4.35). 

Let  zx  =  (24/25,  -7/25)  and  z2  =  (-24/25,  -7/25)  be  two  points  on  the  bound- 
ary 3D,  and  let  P(x,z)  denote  the  Poisson  kernel.  Then 

P(a,  *,)  =  P(b,  z2)  =  lb/221,     P{a,  z2)  =  P(b,  zx)  =  75/29. 

Let  0  =  P(a,  zx)  +  P(a,  z2)  =  P(b,  zx)  +  P{b,  z2),  and  let 

h(x)  =  l(P(x,z1)  +  P(x,z2)). 

Then  h(a)  =  h{b)  =  7.  Further,  P(x0,^i)  =  P(x0,z2)  =  25/51,  so  that 

1 281 8 
h(x0)  =  ^^  *  0.3723456  .  (4.38) 

By  comparing  (4.37)  and  (4.38)  we  see  that  h(x0)  is  strictly  smaller  than 
Rfa(x0)  +  Rfb(x0).  This  proves  that  (Z>,Hinf)  is  not  a  balayage  space. 

Now  we  show  that  ?finf  is  not  a  function  cone  as  defined  in  Section  2.1,  because 
the  axiom  (F3)  does  not  hold. 

Assume,  on  the  contrary,  that  it  holds.  Then  for  every  v  €  Hini,  there  is 
u  £  "Hinf,  such  that  for  every  e  >  0  there  is  a  compact  set  K  in  D  with  v(x)  <  eu(x) 
for  every  x  €  D\  K. 

Let  us  take  v  =  P{-,z),  where  z  —  (1,0)  G  dD.  Then  for  u  as  above  we  take 
c  <  2^oT  ar,d  choose  a  compact  set  K  satisfying  P(x, z)  <  eu(x)  for  every  x  €  D\K. 
Since  u  is  in  HinI,  it  satisfies  for  each  x  E  D  Harnack  inequality 


Then 


»/     x     1-M2       /  x      1-kl2     ,  x     1  1-M2 
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which  implies 


2  <i     for  every  x  G  D  \  K. 


\z-x\2        2' 
We  may  take  x  =  (a,  0)  which  is  m  D  \  K  when  a  is  close  to  1.  Then  the  left-hand 

side  above  is  1,  which  yields  a  contradiction. 

4.5     On  Extremal  Functions 

Let  D  be  the  unit  ball  in  Rn.  Nonnegative  superharmonic  functions  on  D  form 
a  convex  cone.  The  description  of  extremal  functions  in  this  cone  is  well  known: 
The  only  extremal  functions  are  P(-,z),  z  G  dD,  and  Gr>(-,y),  y  G  D,  where  P 
denotes  the  Poisson  kernel  and  Gd  the  Green  function  for  D.  If,  instead  of  the  cone 
of  superharmonic  functions  we  consider  Tiird  on  D,  functions  P(-,z),  z  G  dD,  are 
still  extremal.  On  the  other  hand,  Go(-,y)  are  not  even  in  Ti™1 .  Hence,  they  cannot 
be  extremal.  We  saw  in  Section  2.1  that  sy  =  k(-,y),  y  E  -D,  are  extremal  functions 
in  7iin{.  In  this  section  we  give  some  other  examples  of  extremal  functions.  We  do 
not  know  the  complete  description  of  extremal  elements  in  7iini . 

For  a  point  £  G  D  let  —  £  denote  the  point  on  3D  opposite  to  (,  (0,  ()  =  {\(  : 

0  <  C  <  1}  and  [0,0  =  {A(  :  0  <  (  <  1}.  Recall  the  inequality  (4.1): 

1  -  \x\2 

H^(O)  <  h(x). 

(l  +  |a;|)"    w"    v  } 

Let  z  G  dD;  if  x  G  (0,  -z),  then 

l-|x|2         l-|x|2 


P(x,z)  = 


\z-x\n       (l  +  |x|)n' 
Hence  the  inequality  above  can  be  rewritten  as 

h(x)  >  P(x,  z),    x  G  (0,  -z),  h  G  Ho.  (4.39) 

The  following  lemma  shows  that  unless  h  =  P{-,z),  then  the  inequality  above  is 
strict. 


64 

Lemma  4.2  Let  z  e  dD.  IfheH,  h(0)  =  1  such  that  h  ^  P{-,z),  then 

h(x)  >  P(x,  z),    for  every  x  6  (0,  —z).  (4.40) 

PROOF:  If  h  =  P(-,  £),  where  (  ^  z,  then  for  any  x  £  (0,  —z)  we  have  \x  —  (\  <  \x  —  z\. 
The  formula  for  the  Poisson  kernel  implies  (4.40).  For  an  arbitrary  /iGW,  there  is  a 
probability  measure  n  on  dD  such  that  h(x)  =  fdD  P(x,  ()fi(d().  Using  the  fact  that 
for  some  e  >  0,  n  charges  the  set  {(  €  dD  :  P(x,()  >  P(x,z)  +  e}  and  by  what  was 
shown  above,  it  follows  that  h(x)  >  P(x,z).  □ 

Remark:  With  the  same  assumptions  as  in  the  lemma,  one  can  show  that 
h(x)  <  P(x,z),  for  every  x  E  (0,z). 

LEMMA  4.3  Let  {hn}  be  a  sequence  of  positive  harmonic  functions  in  D  such  that 
hn(0)  >  1  ond  hn(x)  —*  P(x,z)  for  some  x  €  (0,-2).  Then  there  is  a  subsequence 
{hnk}  such  that  hnk  — *  P(-,z)  everywhere. 

PROOF:  The  sequence  {hn}  is  bounded  at  x.  Therefore,  by  Lemma  2.5,  there  is  a 
subsequence  {hnk}  and  a  function  h  (necessarily  harmonic)  such  that  hnk  — »  h  in 
C(D).  Further,  h(0)  >  1  and  h{x)  =  \imkhnk(x)  =  P(x,z).  By  (4.39),  h(x)/h(0)  > 
P{x,z).    Therefore,  h(0)  =  1  and  h(x)  =  P(x,z).    Now  Lemma  4.2  implies  that 

h  =  P(-,z).a 

Lemma  4.4  Let  u  e  H™*,  such  that  u(0)  >  1  and  u(x)  =  P(x,z)  for  some  x  6 
(0,  —  z),  z  6  3D.   Then  u  <  P(-,z)  everywhere. 

PROOF:  Let  u  =  infQ/ia,  ha  G  H.  For  x  as  above,  there  is  a  sequence  {hn}  C  {hQ} 
such  that  u(x)  =  limnhn(x).  Hence  hn(0)  >  1  and  \imnhn(x)  =  P(x,z).  By  Lemma 
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4.3,  there  is  a  subsequence  {hnk}  such  that  limfc  hnk  =  P(-,  z)  everywhere.  Then,  for 
any  y  G  D, 

u(y)  =infft„(y)  <  inf  hnh(y)  <  lim hnk(y)  =  P(y,z). 

Thus  u  <P(-,z).n 

The  next  proposition  shows  that  any  infimum  of  Poisson  kernels  is  an  extremal 
element  in  Hiat. 

PROPOSITION  4.1  Let  {z,  :  i  G  /}  be  any  family  of  points  on  the  boundary  of  D  and 
let 

u  =  inf  P(-.Zi). 

«€/ 

TAen  u  is  an  extremal  element  in  Ti1"1 . 

PROOF:  Assume  that  u  =  f  +  g,  f,g  G  H"**.  Both  /  and  g  satisfy  inequality  (4.1), 
i.e., 

Let  us  fix  i  G  J  and  let  x  G  (0,  — z,-).  Then 

/(*)  +  g(x)  =  u(x)  =  P(x,  *,)  =  i  ~  j|ji  <  /(»)  +  ^(x) 
since  /(0)  +  #(0)  =  1.  Hence, 

/(x)  =  (i  +  M)"/(0)  and  5(:r)  =  (i  +  |x|)^(0)>  x  G  (0'  ~*)-         (4-41) 

Let  /  =  ///(0),  g  =  g/g(0).  Then  /(0)  =  1  and  f(x)  =  P(x,Zi)  for  x  G  (0,-z,)- 
Hence,  by  Lemma  4.4,  /  <  P(;Zi).  Since  this  holds  for  every  i  G  /,  it  follows  that 
/  <  u,  i.e.,  /  <  /(0)u.  Similarly,  g  <  g(0)u.  But,  u  =  f  +  g  <  /(0)tt  +  #(0)u  =  u, 
which  imphes  that  /  =  f(0)u  and  #  =  g(0)u.n 

Let  us  now  fix  a  point  z  G  <9D.  We  show  that  the  function  u  =  P(-,z)  A  6  is 
extremal  for  every  nonnegative  bounded  harmonic  function  h  satisfying  h(0)  =  1. 
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Consider  the  following  sets:  U  =  {x  G  D  :  P(x,z)  <  h(x)},  V  =  {x  g  D  : 
P(x,z)  >  h(x)},  <mdW  =  {x  (E  D  :  P(x,z)  =  h(x)}.  If  there  is  an  x  G  (0,-*) 
(or  x  G  (0,z))  such  that  P{x,z)  =  h(x),  then  by  Lemma  4.2  (i.e.,  Remark  following 
Lemma),  P(-,z)  =  h.  Since  this  case  is  trivial,  we  assume  that  such  an  x  does  not 
exist.  Hence,  (0, -2)  C  U  and  (0,z)  C  V. 

Let  us  assume  that  u  =  f  +  g,  with  f,g  G  Winf.  For  £  G  (0,  —2),  an  argument 
similar  to  the  one  in  the  proof  of  Proposition  4.1  shows  that 

«•> "  (frw/(0)  and  "w  =  (TT|fiF<0)' 

i.e.,  ///(0)  =  P(-,2)  =  g/g(0)  on  (0,-2).  Lemma  4.4  implies  that  ///(0)  <  P(-,z) 
and  g/g(0)  <  P(-,  2)  everywhere  in  D.  Since  /(0)  +  ^(0)  =  1,  this  gives  that  f(x)  = 
/(O)P(x,  2)  and  ^(x)  =  g(0)P(x,  2)  for  x  €  U.  So 

/  =  /(0)«     and    g  =  g{0)u   in  U.  (4.42) 

By  continuity,  this  holds  on  U  C\  D. 

Note  that  on  V  we  have  /  +  g  =  /i.  Hence  both  /  and  g  are  harmonic  on  V. 
Because  of  dt/  D  D  =  dV  D  V  =  W  and  (4.42),  it  follows  that 

f\dVnD  =  f(0)h   and   g\dVnD  =  g{0)h.  (4.43) 

We  would  like  to  show  that  /  =  f(0)h  and  g  =  g(0)h  in  V.  First  we  show 
that  /  and  f{0)h,  considered  as  functions  on  V,  have  same  boundary  values  at  every 
C  G  dV\{z}.  By  (4.43),  this  is  true  on  dV  n  D.  If  C  €  5V  n  5D  and  C  ^  2,  let 
x  — »•  2  in  V.  Since  /i  <  P(-,2)  in  V,  it  follows  that  lim*-^  h(x)  =  0.  But  /  <  h,  so 
lim^,/(x)  =  0.  Therefore,  for  every  C  6  dV  \  {2},  lim^,/(«)  =  lim*_,/(0)A(:r), 
as  a:  stays  in  V.  The  following  lemma  is  needed  to  conclude  that  /  =  f(0)h  in  V.  It 
can  be  easily  derived  using  [11,  p. 118,  (8.4)]. 
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Lemma  4.5  Let  V  be  a  bounded  domain  in  Rn  and  E  a  subset  of  dV  of  harmonic 
measure  zero.  If  u  is  a  bounded  harmonic  function  in  V  such  that  lim^  u(x)  =  0 
for  every  £  €  dV  \  E,  then  u  =  0  in  V . 

We  assume  that  dimension  n  >  2.  Then  the  point  z  is  of  harmonic  measure 
zero.  So  Lemma  4.5  is  applicable  and  it  follows  that  /  =  f(0)h  in  V.  Similarly, 
g  =  g(0)h  in  V.  Hence, 

/  =  /(0)tt     and    g  =  g(0)u   in  V.  (4.44) 

Relations  (4.42)  and  (4.44)  show  that  /  =  /(0)tt  and  g  =  g(0)u  in  D.  So,  the 
following  proposition  is  proved. 

PROPOSITION  4.2  Let  h  be  a  bounded  nonnegative  harmonic  function  in  D  such  that 
h(0)  =  1  and  let  z  6  3D.   Then  u  —  P(-,z)  A  h  is  an  extremal  element  in  /Hiat. 

Remark:  Note  that  the  boundedness  of  h  is  not  essential.  We  can  simply  require 
that  \imx^zh(x)  =  0.  Then  limx_>^/(x)  =  limI_^/(0)/i(x)  for  every  (  €  dV ',  and, 
consequently,  /  =  f(0)h  in  V.  So,  in  this  case  u  is  also  extremal.  Examples  of  such 
harmonic  functions  h  are  functions  of  the  type 


h(x)=  f  P(x,cMdC), 

JdD 


where  fi  is  singular  with  respect  to  the  surface  measure  and  the  support  of  //  is 
bounded  away  from  z. 

Yet,  some  restrictions  on  h  are  needed,  as  the  following  example  shows.  Suppose 
that  hi  is  nonnegative  harmonic  in  D  with  /ii(0)  =  1.  Let  us  define  h  by  h  = 
1/2(P(-,*)  +  hx).  Then  u  =  P(-,z)  A  h  =  \P(-,z)  +  \{P{-,z)  A  hx)t  hence  not 
extremal. 
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Up  to  now,  the  origin  0  was  fixed.  As  mentioned  before,  there  is  nothing  special 
about  0,  so  we  fix  some  other  point  a  6  D,  a  ^  0.  Recall  that  for  z  G  dD, 
La,z  -  {x  e  D  :  sa(x)  =  P(x,z)/P(a,z)}.    So,  if  we  replace  (0,  -z)  by  La<z,  and 
P{-,z)  by  P(-,z)/P(a,z),  all  results  of  this  section  are  still  valid. 

The  result  corresponding  to  Proposition  4.1  is  the  following:  Let  {z,- :  i  G  /}  be 
any  family  of  points  on  the  boundary  dD,  and  let 


u 


=  inf 


iei  P(a,Zi) 

Then  u  is  extremal  in  7iini. 

Let  us  fix  z  G  dD,  and  let  h  be  a  nonnegative  bounded  harmonic  function 
on  D  such  that  h(a)  =  1.  Then  the  result  that  corresponds  to  Proposition  4.2 
says  that  u  =  (P(-,z)/P(a,z))  A  h  is  extremal.  With  this  fact  we  can  dispose  of 
the  condition  that  h(0)  =  1  in  Proposition  4.2.  Indeed,  let  h  be  any  nonnegative 
bounded  harmonic  function  on  D  and  let  u  =  P(-,z)  A  h.  To  avoid  trivialities,  we 
assume  that  there  exists  a  G  D  such  that  P{a,z)  =  h(a).  The  function  h/P(a,z)  is 
1  at  a,  hence  (P(-,z)/P(a,z))  A  (h/P(a,z))  is  extremal  in  Winf.  This  implies  that  u 
is  also  extremal.  Thus  we  can  state  the  following  result. 

PROPOSITION  4.3  Let  z  be  a  point  on  the  boundary  dD  and  let  h  be  a  nonnegative 
harmonic  function  on  D  which  is  bounded  (or  such  that  lim^-..^  h(x)  =  0).  Then  the 
function  u  —  P(-,z)  A  h  is  extremal  in  H,ni . 


CHAPTER  5 
HARMONIC  FUNCTIONS  AND  A  TIME  CHANGE 


5.1     On  Density  of  Kernel  Functions 

In  this  chapter  we  depart  slightly  from  the  general  course  and  investigate  diffu- 
sions on  a  bounded  domain  D  in  Rn  which  have  the  same  class  of  harmonic  functions 
as  the  Brownian  motion  in  D. 

We  start  with  a  simple  but  important  observation  about  the  Poisson  kernel. 

Let  :r0  6  Rn  and  let  0  <  r  <  R.  We  denote  the  Poisson  kernel  for  the  ball 
B{x0,R)byP: 

pM)  =  ir->R2-J*_-f  (5.1) 

where  x  G  B(xq,R),  C,  6  S(xo,R)  =  {z  E  Rn  :  \z  —  x0\  =  R}.  Let  p,  be  a  signed 
measure  on  the  sphere  S(xQ,  R)  such  that 

/  P(x,(Md()  =  0  (5.2) 

JS(xo,R) 

for  every  x  €.  S(x0,r).  The  function  h(x)  —  JSiXom  P(b,  0/^(^0  is  harmonic  in 
B(x0,R)  and  by  (5.2)  it  is  zero  on  S(x0,  r).  The  maximum  and  the  minimum  prin- 
ciples imply  that  h  —  0  in  B(x0,r)  and,  therefore,  h  =  0  everywhere  in  B(x0,R). 
Hence,  the  measure  \i  =  0. 

Let  J  be  the  inversion  through  the  sphere  S(x0,  p)  where  p  =  VrR,  i.e., 

P2 
J(x)  =  x0  +  - — *- — -(x-x0).  (5.3) 

\x  -  x0r 
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The  radius  p  is  chosen  so  that  J  maps  S(x0,  r)  onto  S(x0,  R)  and  conversely.  If  x,  y  £ 
S(x0,r),  then  an  easy  computation  shows  that  \x  —  J(y)\  =  \J(x)  —  y\.  Therefore, 

P{x,J(y))  =  P(y,J{x)).  (5.4) 

Let  v  be  a  signed  measure  on  S(x0,  r)  and  let  \i  =  v  o  J-1  be  its  image  measure 
on  S(x0,R).  Assume  that  for  each  (  (E  S(x0,R) 

(         P(x,()v{dx)  =  0.  (5.5) 

JS(x0,r) 

Then,  for  each  x  £  S(x0,r), 

I  P(x,(MdQ   =     f         P(x,J(y)Hdy)   =    I         P(y,J(x)Wdy)=0 

JS(xo,R)  JS(xo,r)  JS(xo,r) 

where  the  second  equality  follows  from  (5.4),  and  the  last  from  the  assumption  (5.5) 
(since  J(x)  €  S(x0,  R)).  This  shows  that  fi  satisfies  (5.2)  and  therefore  \i  =  0.  So, 
v  =  0,  too.  This  proves  the  following. 

LEMMA  5.1  The  linear  span  of  the  family  {P(-,  C)|5(x0)r)  :  C  €  S(xQ,R)}  is  dense  in 
C(S(x0,r)). 

PROPOSITION  5.1  Let  D  be  a  bounded  domain  in  Rn  and  let  7~((D)  denote  the  family 
of  positive  harmonic  functions  in  D.  If  B(xo,r)  is  a  ball  contained  in  D,  then  the 
linear  span  of  the  family  {h\s(Xo,r)  '■  h  €  li(D)}  is  dense  in  C(S(xo,r)). 

PROOF:  Let  R  >  0  such  that  D  C  B(x0,R)  and  let  P  be  the  Poisson  kernel  for 
B(x0,R).  Each  function  x  t->  P(x,()  is  in  Ti(D).  By  Lemma  5.1  the  linear  span 
of  restrictions  of  such  functions  to  5(x0,  r)  is  dense  in  C(S(x0,  r)).  Hence,  the  linear 
span  of  restrictions  of  functions  in  H(D)  is  also  dense.  D 

Let  8mD  be  the  Martin  boundary  of  D,  and  d°MD  the  set  of  minimal  points  on 
dMD. 
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Corollary  5.1  Let  K(x,()  denote  the  Martin  kernel  for  D,  where  x  €  D  and 
(  G  dMD.  Then  the  linear  span  of  the  family  {K(-,0\s(x0,r)  •  C  £  d°MD)  is  dense  in 
C(S{x0,r)). 

PROOF:  Let  v  be  a  signed  measure  on  S(x0, r)  such  that  for  every  £  €  d°MD, 

I         K{x,()v(dx)  =  0. 

JS(x0,r) 

If  h  is  positive  and  harmonic,  then  there  exists  a  unique  measure  fi  on  the  minimal 
boundary  d%{D  such  that 

h(y)=  f       K(y,CMd()- 
Hence,  by  an  easy  application  of  the  Fubini  theorem,  it  follows  that 

/         h{x)u{dx)  =  0. 

JS(x0  ,r) 

Proposition  5.2  implies  that  v  —  0.  □ 

5.2     Application  to  a  Time  Change 

As  in  the  previous  section,  D  denotes  a  bounded  domain  in  Rn.  Let  X  = 
(Cl^^tjXtiOtfP1)  be  the  Brownian  motion  in  D  killed  while  exiting  D.  After 
being  killed,  the  process  goes  to  "cemetery"  A.  So,  we  consider  the  process  on  the 
extended  state  space  D&  =  D  U  {A}.  The  lifetime  of  X  is  denoted  by  (. 

For  a  Borel  subset  B  of  D&,  let  Tb  denote  the  hitting  time  of  B, 

TB  =  inf{*  >0:XteB}, 
and  Pg  the  hitting  operator, 

PBf(x)  =  Ex[f{XTB):TB<oo}. 
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Let  A  =  (At)  be  a  continuous  additive  functional  of  X  which  is  strictly  increas- 
ing and  finite  on  [0,  (),  and  let  rt  be  the  inverse  of  A: 

rt  =  mi{s  :  As  >  t}.  (5.6) 

Let  X  =  {il,^^T(t)1XT^)j0T(t),Px)-  ^  is  shown  in  [3,  V(2.11)]  that  X  is  a  standard 
process. 

Let  X  =  ($7,  J7,  Tf,  Xt,  0t,  Px)  be  a  standard  Markov  process  on  the  state  space 
D&  with  the  lifetime  £.  We  assume  that  paths  of  X  are  continuous  on  [0,  ().  The 
hitting  time  toBc  D&  and  the  hitting  operator  are  denoted  by  Tb  and  Pb-  To  avoid 
notational  complications,  we  write,  for  example,  Ex[f(XTB)]  instead  of  Ex[f(Xf)]. 

PROPOSITION  5.2  Assume  that  for  every  open  ball  B  =  B(x0,r)  contained  in  D 

PBc(x,-)  =  PBc(x,-)  (5.7) 

for  every  x  €  B.  Then  there  exists  a  continuous  additive  functional  A  of  X  which  is 
strictly  increasing  and  finite  on  [0,£),  such  that  if  r  is  the  inverse  of  A,  the  processes 
X  and  X  (where  X  is  defined  above)  are  equivalent. 

PROOF:  Note  that  the  assumption  (5.7)  means  that  the  exit  distributions  from 
open  balls  are  equal  for  both  processes.  From  the  proof  of  Theorem  1.2  in  [20]  it 
follows  that  the  exit  distributions  from  all  relatively  compact  open  subsets  U  in  D 
are  equal,  i.e.,  P[/c(x,-)  =  PUc(x,-)  for  each  x  6  U.  It  is  easy  to  see  then  that 
P{/c(ar,  •)  =  PVc(x,-)  for  all  x.  From  this  it  follows  that  for  every  compact  subset 
K  of  Da,  Pk{x,-)  =  Pk(x,-)-  The  statement  of  the  proposition  follows  from  the 
Blumenthal-Getoor-McKean  theorem  (cf.  [3,  V(5.1)]).  □ 

Remark:  Whenever  the  conclusion  of  the  proposition  holds,  we  will  say  that  X  is  a 
time  change  of  X . 
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As  in  Chapter  3,  we  say  that  a  function  h  :  D  — *  R  is  harmonic  for  X  if 
Pfcch(x)  =  h(x)  for  every  compact  set  K  in  D,  and  every  x  E  D.  Let  "H  denote  the 
family  of  all  nonnegative  harmonic  functions  for  X.  By  H.  we  denote  nonnegative 
harmonic  functions  in  D.  These  are  precisely  nonnegative  harmonic  functions  for  X. 

PROPOSITION  5.3  Assume  that  T{  —  H.  Then  X  is  a  time  change  of  X . 

PROOF:  According  to  Proposition  5.2,  it  is  enough  to  check  that,  for  any  open  ball 
B  —  B(xQ,r)  contained  in  D  and  any  x  €  B,  the  measures  Pgc(x,  •)  and  Pfic(x,-) 
are  equal.  If  B  denotes  the  closure  of  B,  then  for  h  E.  H,  Psch(x)  =  Pgch(x),  and 
similarly,  since  W  =  H,  pB*h(x)  =  pQCh(x)  for  all  x.  Let  us  fix  x  6  B.  From  the 
previous  observation  it  follows  that  Ps^h(x)  =  pB^h(x)  for  every  h  £  7{.  Since  X 
and  X  are  continuous,  both  measures  Pbc(x,  •)  and  Pbc(x,  •)  are  concentrated  on  the 
boundary  dB  =  S(x0,  r)  of  B.  Hence,  it  follows  from  Proposition  5.1  that  these 
measures  are  equal.  □ 

Let  us  note  that  if  {hn}  is  a  sequence  of  harmonic  functions  for  X  which  con- 
verges in  C(D)  to  a  function  h,  then  h  is  also  harmonic  for  X.  This  was  proved  in 
Proposition  3.1. 

We  assume  now  that  D  is  a  bounded  Lipschitz  domain  in  Rn  and  fix  x0  6  D. 
Let  u>x(A)  denote  the  harmonic  measure  of  a  Borel  subset  A  of  dD  at  x.  Then 
x  h-*  wx(A)  is  harmonic  and  u>x(A)  =  PX(X^  G  A).  For  a  point  z  €  dD,  let  A(z,  r) 
denote  the  "surface  bah"1  on  dD  around  z  (see  Section  3.1  for  a  precise  definition) 
and  let  K  denote  the  kernel  function  for  D  based  at  x0-  Then 

„.       ,       ,.       ux(A(z,2-n))  ,      x 

K(x,z)  =  hm  — v,  A\  '  n    \\,  (5.8 

V  >         n-oou,^o(A(z,2-n))'  V         ^ 

and  convergence  is  uniform  on  compacts.  This  was  proved  in  [15,  Thm.  4.1]  (cf.  also 
[4,  Thm.  3.1]). 
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We  introduce  an  additional  hypothesis  on  X:  We  assume  that  X^_  =  linij^Xt 
exists  in  dD. 

For  a  Borel  subset  A  of  dD,  let 

h(x)  =  PX(X^  6  A). 

If  K  is  a  compact  set  in  D, 

Ex[h{X{TKc))\  =  Ex[Px^c)(Xc.  G  A))  = 

P*[X(TkcHo6kc)-  6  A]  =  P*[*C-  E  A]  =  fe(ar). 
Hence,  a;  i— >  PX(X(-  G  ;!)  is  harmonic  for  X. 

THEOREM  5.1  Let  X  be  the  Brownian  motion  in  a  bounded  Lipschitz  domain  D, 
and  let  X  be  a  continuous  standard  Markov  process  in  D  such  that  X^_  exists  in  dD. 
Assume  that 

PX(X<-  e  A)  =  PX(X{_  G  A)  (5.9) 

for  each  Borel  subset  A  of  dD  and  each  x  G  D.  Then  X  is  a  time  change  of  X. 

PROOF:  Let  z  G  dD  and  x  G  D.  By  (5.9)  and  (5.10), 

„t       v       ,.       m-gA(z,2-"))  P*(XC- e  A(z,2-")) 

Kix.z)  =  hm   „     .    * — — - — —   =    hm  — . 

V'^      n^oop*o(xc_€  A(2,2-"))         ^P^^.eA^^-")) 

Hence,  x  i->  i^(x,  z)  is  harmonic  for  X.  Let  B  =  B(y0,  r)  be  a  ball  in  D  and  y  (z  B. 
Then 

/  K(x,z)PBc(y,dx)  =  K(y,z)=   I         K(x,z)PBc(y,dx). 

JS(yo,r)  JS(yo,r) 

Corollary  5.1  implies  that  harmonic  measures  PBc{y,  •)  and  pBc(y,  •)  are  equal.  Propo- 
sition 5.2  completes  the  proof.  □ 

We  note  that  the  theorem  holds  also  for  NTA  domains.   It  is  enough  to  check 
that  the  equation  (5.9)  is  valid.  But  this  is  proved  in  [16,  Thm.  5.5.]. 
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